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ABSTRACT 

We construct the rest-frame 2-10 keV intrinsic X-ray luminosity function of Active Galactic 
Nuclei (AGNs) from a combination of X-ray surveys from the all-sky Swift BAT survey to the 
Chandra Deep Field-South. We use ~ 3200 AGNs in our analysis, which covers six orders of 
magnitude in flux. The inclusion of XMM and Chandra COSMOS data has allowed us to investi¬ 
gate the detailed behavior of the XLF and evolution. In deriving our XLF, we take into account 
realistic AGN spectrum templates, absorption corrections, and probability density distributions 
in photometric redshift. We present an analytical expression for the overall behavior of the XLF 
in terms of the luminosity-dependent density evolution, smoothed two power-law expressions in 
11 redshift shells, three-segment power-law expression of the number density evolution in four 
luminosity classes, and binned XLF. We observe a sudden flattening of the low luminosity end 
slope of the XLF slope at z > 0.6. Detailed structures of the AGN downsizing have been also 
revealed, where the number density curves have two clear breaks at all luminosity classes above 
log Lx > 43. The two break structure is suggestive of two-phase AGN evolution, consisting of 
major merger triggering and secular processes. 

Subject headings: galaxies: active - galaxies: luminosity function - quasars: general - X-rays: - galaxies 
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1. Introduction 

The AGN/QSO luminosity function and its 
evolution with cosmic time are key observa¬ 
tional quantities for understanding the origin of 
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supermassive black holes (SMBHs) and accre¬ 
tion onto them. This is one of the most im¬ 
portant observational product along with the 
demography of SMBHs, which are made by 
scalin g relations with hos t galaxy properties 
(e.i 


ahn g relations with hos t galaxy properties 
,g. Magorrian et al. 19981: Hdring fe Rfx) 2004 


Gultekin et al . 20091 see TKormendv fc HolbOlil for 


scured (type 1) AGNs except at very high red- 
shifts, imaging surveys at 2 < L[keV] < 10 
with ASCA , XMM-Newton and Chandra enabled 
us to expanded the investigation to include ob¬ 
scured (type 2) AGNs, predominantly within 
the Compton-thin limit (logA^n [cm -2 ] < 24), 
providing a view of the SMBH accretion pro¬ 


recent review) as well as emerging AGN cluster¬ 
ing measurements, which locate typical masses 
of dark matter halos (DMHs) that the AGNs 

cess with much better sampling. Analysis of 
the 2 < i?[keV] < 10 XLF has been made 
by various groups using different samples with 
different approaches and different levels of so- 

reside (e.g., Mullis et al. 2004 Yang et al. 2006t 

Coil et al. 2009| Gilli et al. 

2009| Cannelluti et al. phistication (e.g., Ueda et al. 2003; Cowie et al. 

201C; Krumne et al. 2010, 

20121) or the halo oc- 

2003: La Franca et al. 20051 Silverman et al. 20081; 

cupation distribution (e.g., Mivaii et al 

2011: 

Ebrero et al. 20055" Ycncho et al. 2009^ Aird et al. 

Allevato et al. 2012; Richardson et al. 2013 

Chatteriee et 

21.10; Ueda et al. 20141). Many groups took a sim- 

2013). pie approach and constructed the 2-10 keV XLF 


X ray surveys are practically the most efficient 
means of finding active galactic nuclei (AGNs) 
over a wide range of luminosity and redshift. X- 
ray emission from almost all extragalactic point X- 
ray sources with luminosities of Lx > 10 42 ergs -1 
are considered to originate from SMBH accretion, 
because X-ray emission from the star-formation 
origin (through supernova remnants and X-ray bi¬ 
naries) of galaxies with the highest star formation 
rates typically has a luminosity lower than this 


assuming a simple power-law spectrum without 
absorption corrections. From the combination of 
AGNs from H E AO-1, ASCA and early Chandra 
survey samples, Ueda et al. ( 2003h constructed an 


absorption-corrected 2-10 keV XLF with a full 
./Vh function, as well as a new AGN population 
synthesis model of the X-ray backgro und. This 
work was recently revised by llJeda et ah ( 2014 1 


(IPersic et al .11200-4 iRanalli et al.ll2005t IPtak et al 


2007 1. Enormous efforts have been made by sev¬ 


eral groups to follow up the survey X-ray sources 
with major optical telescopes around the globe, 


(lrerefter, U14), which included both soft and 
hard samples from diverse sources derived from 
ROSAT , MAXI , ASCA , XMM-Newton and Chan¬ 
dra surveys. U14 included sophisticated spectral 
model templates, a revised population synthesis 
model, discussions on possible contributions from 


so that we have fairly complete samples of X-ray 

Compton-thick AGNs as well as the growth of 

selected AGNs ( 

Brandt & Hasingerl 20051 for re- 

SMBHs. A few works have included absorption 


view). One of the most important products of 
such X-ray AGN surveys is the X-ray luminosity 
function (XLF) and its evolution with cosmic time. 
Since X-ray emission marks the accretion activity 
onto the SMBHs, the XLF and its evolution gives 
a bird’s eye view of the moments of SMBH growth 
process. 

The launch of ROSAT and the enormous opti¬ 
cal followup efforts of the X-ray sources detected 
by this satellite enabled us, for the first time, 
to probe the X-ray emitting AGN populations 
to highly cosmological distances. From the mid- 
1990’s to mid-2000’s, the evolution of the XLF 
in the soft band (0.5-2 keV) was investi gated 


2009| h Aird et al. 1 201f)h took a simpler approach 


by numerous authors (e.g., Bovle et al. _ 19931: 


2000a, 2001 


Page et ahl 19971: Ijones et al. 1997t Mivaii et al 


Hasi nger et al. 20051 ). While the 


soft-band surveys found predominantly unob¬ 


in spectral assumptions, but used sophisticated 
statistical techniques and took into account the 
probability density distribution function (PDF) of 
photometric redshifts and incompleteness correc¬ 
tions, with Bayesian and Markov Chain Monte- 
Carlo methods in their model parameter search. 
In an attempt to make the XLF construction as 
free from absorption corrections as possible, S. 
Fotopoulou et al. (2015, submitted to A&A) re¬ 
stricted their analysis to the 5-10 keV band and 
use the full Bayesian approach in obtaining pa¬ 
rameter constraints. 

Most authors fitted their XLFs to simple an¬ 
alytical formulae. As a common feature of most 
such models, the XLF at a given redshift has a 
(smoothed) two power-law form with a shallower 
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slope at the low luminosity end and a steeper 
slope at the high luminosity end. This is quite 
unlike the Scliechter function, which is customar¬ 
ily used to describe the optical luminosity func¬ 
tion of galaxies and has an exponential decline 
at higher luminosities. As for the XLF evolu¬ 
tion, early evolution models with t he pure lumi¬ 


nosity evolution_scheme (PLE; e.g.. iBovle et al 


Il99.il : Page et ahl 19971) have been superseded 
by more complex expressions. Many authors 
have used the Lumino sity-Dependent Density 
Evolution (LDDE) (e.g., SchmjdLfcGreen 1983: 


n. voiuLion [ mur/] i e.g.. locm mcn 66^reen|_[iyoo; 

Mivaii et, al. 112000;' TJeda et, al!2003t Hasinge^etal 


20051 mod el of one form or a nother. lYencho et al 
(2009) and Aird et al. ( 201Clh considered the mod¬ 


els where a smoothed two power-law XLF evolves 
both in the luminosity and density directions, 
keeping the XLF shape the same. These are called 
the (Independent) Luminosity And Density Evolu¬ 
tion (called the ILDE or LADE) models. Because 
the change of at least the low luminosity slope 
with redshift is observed to be present, generally 
the LDDE model is preferred, as critically com¬ 
pared by S. Fotopoulou et al. (2015, submitted to 
A&A). 

When AGN number densities are viewed as a 
function of redshift in different luminosity classes, 


or show “AGN downsizing” (e.g.. l 

Jeda et al. 2003: 

Cowie et al. 2003| Hasinger et al 

2005j) , i.e., the 


peak of the AGN number density appears at lower 
redshifts for low luminosity AGNs than those with 
high luminosities. This was in contrast to an early 
theoretical prediction bas ed on AGN triggerin g 
by hierarchical merging ( Wvithe fc Loebl 2003 1). 


in which it was predicted that lower luminosity 
AGNs peak at higher redshifts (or emerge ear¬ 
lier in the cosmic history). In recent years, there 
have been growing support for the hypothesis that 
the AGN population might be composed of two 
populations that have been triggered by differ¬ 
ent mechanisms. One is a major merger driven 
triggering, which mainly accounts for high lu¬ 
minosity QSOs and dominant in high redshifts. 
The other is a secular process, which might in¬ 
clude fly-by encounters, minor mergers, and/or 
disk instabilitieSjjmd is dominant in intermediate- 


low luminosity Dra per fc Ballantvnel (e.g. l2012h : 
Georgakakis et al.l (e.g. 2014 ). A number of the¬ 
oretical works, involving cosmological N-body 


and/or N-body+hydrodynamical simulations and 
semi-analytical treatments in identifying AGNs, 
have reproduced and/or explained this AGN 
down-sizing e f fect /e.g.. iDi Matteo et al. 2003: 
Hojokins^^aL AlOS^JCroton 2009: Marulli et al 


2008: Degraf et alJI^OloT Fanidakis et al. 20121). 

Some models also include “radio m o de” ( a.k.a. 
“hot halo” mode). iFanidakis et al. (2 012 ) put 
emphasis on dividing the AGN activity into “star- 
burst mode”, which includes SMBH accretion fol¬ 
lowing both merger-driven and disk instability 
driven starbursts, and the “hot halo” mode. By 
also including empirical obscuration effect, they 
reproduced the AGN down sizing ef fe ct. W ith 
this general scheme, Fanidakis et al. ( 20131 ) ex¬ 
plained the tendency that intermediate X-ray lu¬ 
minosity AGNs are associated with more mas¬ 
sive da rk matter halos (DMHs) than luminous 


QSOs (Mivaii et, al. 2007 : Krumne et al. 201f)t 


Allevato et al. 20111 Krumpe et al.l 20 12L In or 


cler to critically compare these lines of theories 
and observation, refinements in calculations of the 
XLF and investigating its detailed evolutionary 
behavior continues to be an important observa¬ 
tional task. 


In this work, we use AGNs selected at E > 2 
keV from a collection of surveys. Our aim is to 
construct and investigate the detailed behavior 
of the XLF. In this particular work, we make 
our best effort to construct the luminosity func¬ 
tion of unabsorbed and absorbed AGNs within 
the Compton-thin limit, i.e. absorbing column 
densities of logiVnlcm -2 ] < 24. Above this col¬ 
umn density, Compton scattering causes photons 
to travel along longer paths within the absorbing 
medium and are subject to much higher chance 
of photo-electric absorption. In such Compton- 
thick (CTK) AGNs, X-ray emission at E < 10 
keV, which is observable with ASCA , XMM- 
Newton , and Chandra, is highly suppressed un¬ 
less they are at very high redshifts. Thus AGNs 
detected in the E < 10 keV surveys are consid¬ 
ered to be dominated by Compton-thin (CTN) 
AGNs, with the exception of a small number of 
X-ray sources, for which X-ray emission is domi- 
nated by in scattered and reflected components 


(Brunner et al. 

2008: Brightman & Ueda 2012: 

Brightman et al. 

2014). In the local universe. 


higher energy X-ray surveys (20 < E < 200 keV) 
such as those available from Swift BAT (e.g., 
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Tueller et al. 

2010: Aiello et al. 

2012}), INTE- 

GRAL (e.g., 

Beckmann ct al. 2009 

: Treister et al. 


AGNs 24 < log Ah [cm 2 ] < 25 and quantify their 
space density and em issivities. NuSTAR (e.g., 


Alexander et al .1 120131) is extending this to higher 


redshifts (z ~ 0.5 — 1). 

We present here a new estimate of the X-ray lu¬ 
minosity function evaluated at the traditional 2- 
10 keV rest frame. Complementary to a recent 
work by Ueda et al. ( 201411 . whose emphasis is on 


the X-ray population synthesis of the X-ray back¬ 
ground as well as discussions on the role of the 
CTK AGNs, our emphasis is on the detailed be¬ 
havior of the shape and evolution of the XLF as 
well as discussions of possible systematic errors. 
In particular, in order to assess the uncertainties 
due to photometric redshifts, we include the PDFs 
of photometric redshifts (photo-z’s) into analysis 
and compare with the results obtained by using 
single best-fit photo-z’s. We focus on the X-ray 
luminosity function of unabsorbed and Compton- 
thin (CTN) AGNs because of the wide availability 
of E < 10 [keV] AGNs. Also X-ray luminosity Lx, 
which is the independent variable of the XLF pre¬ 
sented in this paper, is primarily defined as the 
intrinsic 2-10 keV luminosity, i.e., before absorp¬ 
tion and not including the reflection component. 


In this work, the sample is greatly expanded 
by the add ition of the datas ets from XMM- 
COSMOS (Hasineer et al.l2007ll and C-COSMOS 
( Elvis et al. 20091) . which are parts of the COSMi c 
evolution Survey (COSMOS) (IScoville et al.ll2007l) . 


In Sect. [21 we explain our COSMOS sample, in¬ 
cluding the construction of our combination of the 
XMM -COSMOS and C-COSMOS datasets. The 
details of the selection criteria for other samples 
from the literature are explained in Sect. [3] At 
the end of this section, we explain our incomplete¬ 
ness correction and derive the A(> S) function 
for extragalactic X-ray sources in each sample. 

In sect. [H we explain the computation of our 
estimated XLFs. Model parameter estimations us¬ 
ing a Maximum-likelihood fitting method, which 
also include absorption and the probability den¬ 
sity distributions (pz-PDFs) of photometric red¬ 
shifts, are explained. This section also covers the 
-Nobs/Amdi estimation of the binned XLF in the 
presence of the absorption and pz-PDFs. 

Sect. [5] presents our global expression of 


the Luminosity-Dependent Density Evolution 
(LDDE). In Sect. EH we take a closer look at 
detailed behavior of the XLF in redshift divided 
shells and the number density/emissivity evolu¬ 
tion in divided in luminosity classes. We also 
present binned XLFs. The results are discussed 
in Sect. [7] and a concluding summary is given in 
Sect. 0 

Throughout this work we use a Hubble con¬ 
stant H 0 = 70/i7o kms _1 Mpc _1 n m = 0.3, and 
= 0.7. The hn dependence is explicitly stated. 
If units are omitted, X-ray luminosities (L with 
any sub/super scripts) are measured in units of 
/i^jfergs -1 , and column densities ./Vh in cm -2 . 
The symbol log signifies a base-10 logarithm and 
In a natural logarithm. 


2. The COSMOS Sample 


2.1. The XMM- COSMOS Sample 

As a part of the COSMOS collaboration, 
the XMM -COSMOS program (PI—0. Hasinger 
Hasineer et al.1 12007 ) included observations the 
entire 2.14 deg 2 COSMOS field with the EPIC 
camera on board the XMM-Newton observatory 
over three observing cycles, AO-3, 4 and 6. X-ray 
point source catalogs have been produced using 
a sophisticated source detection and charact eriza- 
tion procedure (jCappelluti et all2007l [20091 '). The 
procedure was developed as a part of the XMM- 
SAS packageQ- In the combined dataset, which is 
composed of 55 XMM-Newton observations with 
a total exposure of ~1.5 Ms, 1887 unique X-ray 
sources have been detected. Optical counterparts 
and their multiwavelength properties were pub¬ 
lished for the X-ray source det ected in the first 12 
XMM-Newton observations bv lBrusa et al.l ( 20071) 
and all 53 s u cces sful observations made before 
2007 by Brusa et al. ( 2010l ) respectively. 

In this work, we use X-ray sources detected in 
the hard (EPIC energy channel between 2 and 8 
keV) band in the 53 XMM-COSMOS observation 
data with a maximum-likelihood of ML > 10, 
which corresponds to the f alse detect ion probabil¬ 
ity of e -10 ~ 4.5 x 10 -5 . iBrusa et~ ( 2010 ) pub¬ 
lished a multiwavelength catalog for the 53-field 
XMM-COSMOS sources with a primary 2-10 keV 
flux cut of Sx > 3 x 10 -15 ergs -1 cm -2 , which 


J http: //xirnn. esac. esa. int/sas/ 
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is our XMM- COSMOS base catalog. After this 
flux cut, we are left with 923 hard-band selected 
XMM- COSMOS sources. We use the EPIC 2-8 
keV countrate to 2-10 keV flux energy conversion 
factor (hereafter referred to as ECF) obtained by 
assuming a power-law spectrum with photon index 
r = 1.7 absorbed by a Galactic absorbing co lumn 
density of IVh = 2.5 x 10 2O cm -2 ( Cappelluti et al 
2009). 


2.2. The Chandra-COSMOS Sample 

With its unprecedented spatial resolution, the 
Chandra X-ray observatory (CXO) provides su¬ 
perb point source sensitivity in deep surveys. 
The Chandr a COSMOS prog ram (C-COSMOS, 
PI=M.Elvis; Elvis et al. 2009il is a large Chandra 
program with a total exposure of ~ 1.8 Ms with 
49 observations with the ACIS-I instrument. The 
observed fields are arranged in 7x7 overlapping 
tiles, covering a central 0.9 deg 2 of the COSMOS 
field. Inclusion of the information obtained by the 
recent Chandra Legacy COSMOS survey (2.5 Ms, 
PI=F. Civano) will be a topic of a future work. 
A total of 1781 unique X-ray sources have been 
detected with a likelihood threshold of ML > 10.8 
in any of the full (0.5-7keV), soft (0.5-2 keV) or 
hard (2-7 keV) bands, out of which 1017 are hard 
band detected ( Elvis et all 120091 : IPuccetti et al. 
2009lh _ 

In Elvis et al. ( 2009h . the 2 10 keV flux is 


given for an ECF calculated assuming a T = 1.4 
power-law. To treat the sample in the same way 
as XMM-COSMOS, we have converted the 2-10 
keV flux using the ECF assuming T = 1.7 for 
our further analysis. Unlike the case of XMM- 


COSMOS, Civano et al. ( 2012l f did not 


impose 

any further flux limit in their multiwavelength 
identification catalog. However, we have imposed 
a 2-10 keV flux (r = 1.7) limit of Sx > 1-6 x 
10~ 15 ergs -1 cm~ 2 , below which the log A — log S' 
curve spuriously rises due to the Eddington bias. 
This removes the 24 faintest hard-band detected 
so urces. Note that t he log A—log S curve in Fig. 9 
of Elvis et al. ( 2009 1 effectively plots down to this 
flux level. This leaves 993 C-COSMOS hard-band 


sources. 


2.3. The Combined XMM-Chandra COS¬ 
MOS Sample 

2.3.1. Sensitivity map combination 


In order to take advantage of both the XMM 
and C-COSMOS surveys, we have made a com¬ 
bined sample, which we call the XC-COSMOS 
sample, by selecting X-ray sources from the more 
sensitive survey out of the two at any given source 
position. For this purpose, we have compared 


the hard X-rav sensitivity map for the 

53-field 

XMM-COSMOS survey ( 

Cappelluti et al 

2009b 

with that for C-COSMOS (Puccetti et al. 

2009). 


The XMM-COSMOS sensitivity map has been 
truncated at the lowest limiting flux of the opti¬ 
cal catalog described above. In order to compare 
both maps on the same grounds, we scaled the C- 
COSMOS sensitivity map using the ECF for the 
photon index of T = 1.7. We further smoothed 
the C-COSMOS sensitivity map with a Gaussian 
filter with cr = 16" and combined with the XMM- 
COSMOS sensitivity map by choosing the smaller 
limiting flux at each point. We have also gener¬ 
ated a mask file that indicates whether the XMM- 
COSMOS or C-COSMOS source should used at 
each position. In effect, the C-COSMOS sources 
should be used in most of the region covered by 
C-COSMOS, except for small areas near corners 
and edges. The combined sensitivity map is shown 
in Fig. m In the combined sample, there are 
1400 extragalactic point sources (mostly AGNs), 
18 Stars, and 29 unidentified sources. Overall the 
fraction of sources identified with AGNs/galaxies 
with spectroscopic or photometric redshifts and 
stars is 98%. Note that the XC-COSMOS AGNs 
comprises as many as ~ 43% of all AGNs in our 
global AGN sample defined below. 


One of the greatest advantage of the COSMOS 
survey is the availability of multi-epoch, homog¬ 
enized, deep 31 band photometry, covering from 
the UV to the MIR wavelengths, including 12 in¬ 
termediate bands from Subaru (ITaniguchi et al. 
2007). These data, along with the correction 


for variability, and the use of hybrid AGN tem¬ 
plates, have allowed the COSMOS survey to 
reach a mean photometric redshift accuracy of 
^Az/fi+zspec) ~ 0-015 with a low fraction of out¬ 
liers (~5%) for b oth the XMM and Cha ndra 
COSMOS samples ( Salvato et al. 20091 201 lib In 
this work, the newest photometric redshifts by 
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Salvato et aL (2011) have been used when there is 


no available spectroscopic redshift. The associated 
probability density distributions of photometric 
redshifts are also used. In addition to the spec¬ 
troscopic redslrifts included in Brusa et~ahl ( 201Clh ■ 
newly available proprietary spectroscopic redshifts 
from spectroscopic survey programs within the 
COSMOS consortium have been included. 


XMM- & C-COSMOS Combined Sensitivity Map (2-10 keV) 
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Fig. 2.— The sample used in our analysis in the 
redshift-Observed X-ray luminosity space. The 
surveys (or a group of similar surveys) are shown 
in different colors/symbols as labeled. (Color ver¬ 
sion of this figure is available in electronic version 
only.) 
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Fig. 1. The XMM-COSMOS/C-COSMOS com¬ 
bined sensitivity map. The scale indicated by a 
solid line corresponds to 1 degree on the sky. The 
numbers at the tick marks in the gray scale bar 
at the bottom of the figure indicate the base-10 
logarithm of the limiting 2-10 keV flux in units of 
ergs~ x cm~ 2 assuming af = 1.7 power-law. 


3. The Global AGN Sample 
3.1. The Swift BAT Survey Sample 


are evaluating intrinsic XLF at 2-10 keV, we need 
spectral template for unabsorbed AGNs to con¬ 
vert from 15-55 keV to 2-10 keV. In our simplest 
treatment, we assume a power-law spectrum with 
a photon index of T = 1.7, which is approximately 
(within ~ 4%) equal to the corresponding ratio 
for the unabsorbed AGN template spectrum (see 
Sect. 14.311 . Under this assumption, we can convert 
the 15-55 keV flux to 2-10 keV flux by multiply¬ 
ing by a factor of 0.705. I n the limiting flux-su rvey 
area curve from Fig. 1 of Aiello et al. ( 20121) . the 
flux is converted to the 2-10 keV value in this 
way. In our full treatment in Sect. 14.31 we use 
a luminosity-dependent effective photon index for 
the conversion. 


Based on the Swift B AT slew scan data accu¬ 
mulated over 60 months, Aiello et al. (2012) pro¬ 
duced a catalog of AGNs detected in the 15-55 
keV band, providing the current most complete 
census of nearby bright AGNs over the entire sky. 
These include unabsorbed and absorbed AGNs up 
to moderate CTK AGNs. For CTN AGNs, photo¬ 
electric absorption affects very little to the 15-55 
keV flux and thus there is no need for compli¬ 
cated absorption corrections. However, since we 


3.2. ASCA LSS/MSS 

The Advanced Satellite for Cosmology and As¬ 
trophysics (ASCA) was the first X-ray observatory 
with focusing optics capable of making spectro- 
copic imaging in the 2-10 keV band. The com¬ 
bination of ASCA Large Sky Survey (ALSS) and 
the ASCA Medium Sensitivity survey (AMSS) fills 
the gap between all-sky surveys and deeper sur¬ 
veys with XMM-Newton and Chandra, occupying 
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a unique regime in the z — Lx space at E > 2 keV. 

The ALSS covers a continuous area of 5.5 deg 2 
near the North Galactic Pole with a flux limit 
of Sx ~ l x 10~ 13 |’erja; s~ 1 cm~ 2 1 (2 - 10 keV) 
( Ueda et al. 1998 . 19991) . Thirty AGNs detected 


wit h the SIS instrument a re completely identified 


by ( Akivama et al. 2000l) . The AMSS is based 


on a serend ipitous X-ray survey with the GIS 
instrument! Ueda et al.ll200~l . 2005 ) and has a sur¬ 
vey area of ~ 85 deg 2 . The flux limit ranges from 
5 x 10 -12 to 3 x 10“ 13 ergs -1 cm -2 . The identifi- 
cation catalogs in the northern part of the AMSS 
( Akivama et al .1120031 AMSSn;) and the southern 
part (AMSSs; Akiyama et al. in prep.) contain 
74 and 20 spectroscopically identified non-blazer 
AGNs respectively. Three X-ray sources are left 
unidentified and the identification completeness is 
97% for the combined ALSS and AMSS sample. 

The 2-10 keV fluxes of these sources have been 
obtained from the ASCA GIS count rate assum¬ 
ing a T = 1.7 power-law using the GIS response 
matrix. 

3.3. The XMM-Newton Hard Bright Serendip¬ 
itous Sample (HBSS) 

The Hard Bright Se rendipitous Sample (HBSS 
Della Geca et al.l 12004 ) is a subsample detected 
in the 4.5 - 7.5 keV band from the larger sur- 
vey of the XMM -Newton Bright Survey (XBS). 
Della Coca et al. ( 2008 ) defined a complete flux- 
limited sample of 67 sources with MOS2 count 
rates larger than 0.002 [ctss -1 ] (4.5 - 7.5 keV) 
over 25 deg 2 of the sky. The MOS2 count 
rates in the HBSS were converted to observed 
fluxes in the 2-10 keV band, assuming a pho¬ 
ton power-law index of 1.7 (using an ECF of 

cts -1 . The correspond- 
erg s _1 cnT 

sources are completely identified spectroscopically 
except for two, making the identification com¬ 
pleteness of ~ 97%. The 65 identified sources 
consist of 62 extragalactic point sources (mostly 
AGNs), 2 stars and 1 cluster of galaxies. 

3.4. The XMM-Newton Medium-sensitivity 
Survey (XMS) 

Th e XMM-Newton Med ium-sensitivity Survey 
(XMS lBarcons et ah 2007 ) is a serendipitous X- 
ray survey and optical identification program of 


1.044 x 10 11 erg cm 2 

ing flux limit is 2.1 x 10~ 13 erg s _1 cm~ 2 . These 


sources with intermediate X-ray fluxes discovered 
in 25 XMM-Newton high Galactic latitude fields 
covering a sky area of ^3 deg 2 . For this analysis 
the XM S-H sample, sel e cted i n the 2-10 keV band, 
is used. iBarcons et al. ( 20071) assumed a spectral 
index of T = 1.7 for the count-rate to flux con¬ 
version for the hard band, and therefore we use 
the 2-10 keV flux as shown in their dataset. In 
our analysis, we use the the sub-sample of XMS-H 
defined bv lHasinger ( 2008 ). 

The original Barcons et al. ( 20071) sample con¬ 
tains 159 sources with 2-10 keV fluxes > 3.3 x 
10~ 14 [ergcm~ 2 s" 1 ] with a spectroscopic identifi¬ 
cation fraction of 83% (27 unidentified sources). 
However, the actual spectroscopic comp leteness 
limit varies from field to field. Therefore HasingeJ 
(2008 ) selected a subsample that is comprised of 
further flux cuts on a field-by-field basis. This 
way a cleaner XMS sample could be defined, in¬ 
cluding 128 extragalactic point X-ray sources and 
almost half the number of unidentified sources, i.e. 
achieving an identification fraction of 91%. The 
survey solid angle_has been corrected accordingly. 
See Sect. 2.4 of Hasinger ( 2008 ) for further details. 


3.5. Hellas2XMM 


HELLAS2XMM ( Baldi et al.ll2002 ) is a serendip¬ 
itous survey based on suitable XMM-Newton 
pointings (complementary to the XMS). The as¬ 
sumed spectral in dex fo r the ECF for the 2-10 
keV flux in iBaldi et al.l ( 20021) was T = 1.7, and 
therefore we use their flux as provided. 

As in the case of the XMS, we use the sub¬ 


sample compiled bv lHasingeii (|2008l ). iFiore et al 


(12003!) presented optical identifications and spec¬ 
troscopic redshifts for 122 sources selected in the 
2-10 keV band in five XMM-Newton fields, cover¬ 


i ng a s urvey area of 0.9 deg 2 . Later iMaiolino et al 


(2006) presented additional redshifts from VLT 
NIR spectroscopy of optically extremely fai nt ob¬ 


j ects f or two ad ditional sources . BothlFiore et al 


( 2003 ) and later Mignoli et al.l ( 2004 ) tried to es¬ 
timate redshifts for th e remaining, optically faint, 

pub¬ 


unidentified sources. Cocclria et al 


r ticaiiy 
( 20071 ) 


lished photometry and spectroscopic redshifts in 
five additional HELLAS2XMM fields providing 
59 new redshift identifications for the sample of 
110 new sources with an additional solid angle of 
0.48deg 2 . As explained in Hasinger ( 2008 ). further 
flux cuts were made on the field-by-filed basis. In 


7 























































































addition, different flux cuts were applied to inner 
and outer regions of each field. The redshift com¬ 
pleteness in this subsample of HELLAS2XMM is 
thus 84%. In this paper, we furth er raise the flux 
limit from l.OxlO -14 (adopted bv lHasinger 2008 1 


to 1.63x10 14 ergs 'em 2 to make the identifica¬ 
tion 90% complete, leaving 115 AGNs. 

3.6. Lockman Hole (XMM) 


Because of the extremely low Galactic column 
density (Ah = 5.7 x 10 19 cm -2 ), Lockman hole 
has been selected as a target of X-ray deep survey 
observations, includi n g a XMM-Newton survey 


( Hasinger et al. 20011 Brunner et al.1 2008 b The 


base X-ray so u rce ca talog for this work is from 


Brunner et al. ( 20081) from 18 individual point¬ 


ings, twelve of which are centered near the same 
pointing direction and the remaining pointings 
are sp read out over a bout ~ 30' in right ascen- 
Brunner et al. ( 20081) concentrated on the 
central circular area with a radius of 15' 


sion. 

deep 


(= 0.20deg 2 ). Their catalog contains 409 X-ray 
sources, out of which 266 are detected in the hard 
(2-4.5 keV) band at the maximum-likelihood (ML) 
value of 6.0 or above. Based on simulations, they 
showed that the sources above this ML threshold, 
in combination with the flux-sensitivity curve in 
their Fig. 5, reproduces the log A — log S accu¬ 
rately. 

Their listed 2-10 flux is converted from the 
EPIC 2-4.5 keV count rate assuming a T = 2.0 
power-law spectrum. For our purposes, we have 
converted their T = 2.0 based 2-10 keV flux to the 
T = 1.7 based flux using a combined response ma¬ 
trix of the XMM-Newton EPIC PN and two MOS 
detectors. This corresponds to a multiplication of 
the r = 2.0 based flux by 1 .16. The li miting -flux 
- area curve in Fig. 5 of Brunner et alJ ( 20081) has 


also been converted to the T = 1.7-based flux. 

Of the 266 hard X-ray selected sources, 60% 
have published or unpublished spectroscopic red- 
shifts (including 5 stars). Inclu ding p hotomet- 
ric redshifts bv IFotopoulou et al. ( 20121 ). 98% are 
identified with either spectroscopic or photomet¬ 
ric redshifts. Only 4 are left unidentified. This re¬ 
sults in 257 hard-band selected extragalactic X-ray 
point sources. Probability density distributions of 
the photometric redshifts are also available, which 
we use in Sect. 14.21 


3.7. CLASXS and CLANS 

The catalogs from two intermediate-depth 
Chandra surveys near the Lockman Hole, the 
Chandra Large Area Syno ptic X-ray Survey 
(CLASXS; lYang et al 20041) and t he Chandra 
Lock man Area North Survey (CLANS lTrouille et al.1 
120081) are also included in our analysis. The X- 
ray source identifications, and limiting-flux - sur¬ 
vey area curves of both surveys (as well as those 
of the C DF-N; see bel o w) ar e conveniently in¬ 
cluded in iTrouille et al.1 (20081) as a part of their 
“Opt-X” project. Additional spectroscopic red¬ 
shifts of CLANS sources can also be found in 
( 20091 ). Their tables contain 2-8 


Trouille et al 


keV fluxes based on the 2-8 Chandra ACIS-I count 
rates with the ECFs calculated with variable T 
and ACIS-I response matrices based on old cali¬ 
brations. Therefore we have converted from their 
2-8 keV count rate to 2-10 keV flux assuming 
r = 1.7 and with ancillary response files (ARFs) 
we have re-created for the ACIS-I aim point using 
the Chandra CALDB 4.5.7 (the latest calibration 
as of our analysis). For each of the CLASXS and 
CLANS surveys, the applied response is the expo¬ 
sure weighted average of the responses at the aim 
points of the observations belonging to the survey. 
The flux i n their limiting flux -survey area curves 
(Fig. 5 in ITrouille et al.ll2008l ) have been adjusted 
accordingly. The curves for the probability of de¬ 
tection of 30% are used, because they reproduce 
the lo g A — logS relation in the simulations by 
lYang et al ( 2004 ). Because their curves only ap¬ 
ply to the sources with offaxis angles of < 8', we 
use sources that meet this criterion. Also, we only 
include X-ray sources that are detected at a signal- 
to-noise ratio of > 3 in the hard band in order to 
be compatible with their flux-survey area curve. 
Moreover, in order to be compatible with there 
flux-limit - area curve, they also imposed offaxis- 
dependent flux limits to exclude some sources for 
their log A — log S' calculations. These imposed 
limits cannot be reproduced with the published 
information. However, this effect does not signifi¬ 
cantly affect our analysis as described below. 

In order to minimize the effects of the “com¬ 
pleteness correction” (Sect. 13.121) . we impose at 
least 90% identification completeness for each sur¬ 
vey. In order to achieve this, we have excluded the 
sources with 2-10 keV flux (T = 1.7 based) of Sx < 
6.5 x 10~ 15 and S < 6.5 x 10 -15 ergs -1 cm -2 , leav- 






































































ing 250 and 119 AGNs for CLASXS and CLANS 
respectively. At these flux limits, the survey area 
is already ~ 95% of the total geometric area with 
offaxis angle < 8'. Thus the effects of offaxis- 
dependent flux limits are minimal. 


3.8. Chandra Deep Field North (CDF N) 


Our base X-r ay source catalo g for the 2Ms 


CDF-N is from Alexander et aJj (20031 ). As in 
the cases of CLASXS and CLANS, we have used 
the hard-band detected sources that have signal- 
to-noise ratio of 3 or better to match with the 
flux-limit survey area curve shown in their Fig. 
19. In order to convert their 2-8 keV count rate 
to 2-10 keV flux, we have re-generated the ancil¬ 
lary response file for an exposure-weighted aver¬ 
age o^_all_the_CDF-N_Chandra observations used 
by Alexander et al- (20031 ). assuming a h = 1.7 
power-law and using CALDB 4.5.7. Since their 
flux-limit survey area curve is based on T = 1.4 
and old calibration, we have adjusted their curve 
for our flux definition. 

We use the opt ical identifications listed in 


Trouille et al. ( 20081) . which include redshifts from 


the literature as well as from their own program. 
We have further imposed a lower flux cut of 
Sx > 8.9 x 10 -16 ergs -1 cm -2 to ensure an identi¬ 
fication completeness of at least 90%, leaving 182 
extragalactic X-ray point sources. 

3.9. Chandra Deep Field South (CDF—S) 

An extensive X-ray source and identification 
catalog of the 4 Ms Chandr a Deep Field South 
was published by IXue et all ( 201ll) . which also 
includes a list of previously published spectro- 
scopic and photometric redshifts. Subsequently 
Lehrner et al. ( 20121 ) generated limiting flux 
solid angle curves and the log N — log S relation. 
This work took into account flux probability dis¬ 
tributions and recovery functions for those sources 
with a threshold probability of false dete ction of 
Pthresh = 0.004, as adopted bv lXue et al. ( 2011 ). 

We have chosen hard X-ray sources that have 
been detected in the 2-8 keV band from their cat¬ 
alog. Their count-rate to flux conversion is based 
on the effective photon index T from the ratio of 
the hard and soft band count rates, and therefore 
the ECF is variable. For each source, we have con¬ 
verted their 2-8 keV flux (based on their estimated 


T) to 2-10 keV flux based on T = 1.7 using the aim 
point response matrix of one of the CDF-S ob¬ 
servations (OBSID=12049). As shown below, the 
value of the count-rate to flux conversion factor is 
subject to the variation of instrumental response 
over time as well as calibration updates, but the 
relative change of conversion factors among differ¬ 
ent spectra is insensitive to ACIS-I response ma¬ 
trix used. The 2-10 keV flux (after the conversion) 
of the faintest object amongst the hard-band de¬ 
tected sources was 8.3xl0 -17 ergs -1 cm -2 . 


We have used the 2-8 keV band limiting flux 
- soli d angle curve in Fig. 2 of iLehmer et al. 
(120121) . Their curves are drawn for T = 1.0 and 
we have converted their relation to T = 1.7 using 
the method described above. 


We have further included new sp ectroscop ic 


redshifts that are not included in the Xueet al 


( 201ll ) catalog. In the catalog of Xue et al. ( 2011 ). 
there are two objects with photometric redshifts 
above z n h > 7. The secon d solution of these ob¬ 
jects by ILuo et ahl ( 2010) are both at z p h ~ 3. 
However, a new spectrocopic redshift obtained for 
one of them (XID=28) is z sp = 0.69. Thus it is 
possible that the other also suffer a similar catas¬ 
trophic photometric redshift error also. Thus we 
have excluded the remaining z p h > 7 object from 
our sample. The overall identification complete¬ 
ness is 97% (including photometric redshifts) and 
the spectrocopy completeness is 71%. 


3.10. Removing Compton-thick AGNs 

The 2-10 keV selected samples, especially at 
low redshifts, are highly selected against CTK 
AGNs. Even the Swift BAT sample selected at 
15-55 keV suffers from suppression of the CTK 
population, because of the Compton scattering it¬ 
self which results in subsequent lowering of the 
photon energy as well as longer paths in absorb¬ 
ing medium. However, there always is some spill 
over of the CTK population to the sample defined 
here. There are two approaches in treating the 
situation. One is to include the CTK population 
into the model and make model fits to the sam¬ 
ple that may include CTK AGNs. In other words, 
the amount of CTK AGN spill-over is estimated 
by the model. The model of U14 introduced the 
parameter /ctKj which is the ratio of the number 
density of the CTK AGNs (24 < log 1V H < 26) to 
that of absorbed CTN AGNs (22 < log Ah < 24), 
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and they then assumed the same /ctk at all lu¬ 
minosities and redshifts. They showed that this is 
at least consistent with various observational es¬ 
timates of the CTK AGN population in the lit¬ 
erature. However, U14 also estimated that the 
observational constraint on this parameter was 
0.5 < /ctk <1.6, carrying a factor of three un- 
certainty. This is con sistent with recent results by 


Buchner et al.l (1201 5). An alternative approach is 


to remove known CTK AGNs from the sample on 
a best effort basis and construct an XLF model 
that represents only the CTN AGNs. We take the 
latter approach. 

There are a number of studies aimed at identi¬ 
fying CTK AGNs among the AGN samples used 
in this paper. Because the Swift BAT sample is se¬ 
lected in the 15-55 keV band, it has some sensitiv- 


i ty to CTK AGNs with log A/j < 25. lAiello et al 


( 20121 ) list 18 known CTK AGNs from their pre¬ 
ceding work of Burlon et al.l ( 2011 ) and other lit- 
erat ure. These have bee n excluded from our anal¬ 
ysis. lAiello et~ah ( 20121) estimated that a few un¬ 
known CTK AGNs might still be present in this 
sample, the effect of which on our analysis is negli¬ 
gible. In their XMM-Newton Lockman hole sam¬ 
ple, 11 AGNs meet the X- ray color-color criteri a 


of CTK AGNs defined bv I Brunner et al.l ([2008). 


None of them is in our hard -band selected s ample 


above our luminosity cut. ILanzuisi et al.1 (12015 ) 


selected 10 CTK AGNs based on spectral fits for 
the XMM-COSMOS sources that have 30 counts 
or more in the full 0.3-10 keV band. None of 
them are in the XMM part of the XC-COSMOS 
AGN sample. Note that all but one XMM-Newton 
AGNs in the XMM part of the XC-COSMOS sam- 
ple have 0.5-10 keV counts greater than 30, thus, 


Lanzuisi et al. ( 20151) would have identified almost 


all CTK AGNs in it if any were present. One AGN 
in our CLANS sample i s among the fi v e X-ra y se¬ 
lected CTK AGNs by Polletta et, al.l ( 2007 ) and 
this has been excluded for our analysis. Nine of 
the CDF-N AGNs are among the CTK AGN can- 
didat es d erived from the X-ray spectral analysis by 
Georgantopoulos et ali ( 2009i) . These are also ex¬ 


cluded from our analysis. Georgantopoulos et al.l 
( 2009 ) performed spectral analysis of the CDF- 
N AGNs with a 2-10 keV flux of > 1 x 10 -15 


[ergs -1 cm -2 ], which is approximately the same 
as our flux cut to the CDF-N so urces to meet our 
completeness criteria. Recently iBriglitman et al 


(|2014P published a table of extensive X-ray spec¬ 
tral analysis results of the X-ray sources in the 
CDF-S, AEGIS-XD and C-COSMOS surveys with 
an intention to identify CTK AGNs. Because of 
the low background level of Chandra surveys, they 
attempted spectral analysis to the sources down 
to 10 C handra ACIS counts. F ollowing the cri¬ 
teria by Briglitman et al.l (2014), we consider as 
probable CTK AGNs those that satisfy the both 
of the conditions; the best-fit value of logA/j is 
greater than > 24 and the lower bound of its 
90% confidence error range is > 23.5. We ex¬ 
clude these from from our analysis also. The num¬ 
ber of excluded sources by the CTK AGN criteria 
are 24 for the C-COSMOS and 23 for the CDF-S 
samples respectively. The full band counts of al¬ 
most all faintest sources above our flux cuts (see 
above) well exceed 10 counts for both C-COSMOS 
and CDF-S samples. Thus CTK AGNs that may 
be still remain in our sample can be neglected. 
The CTK AGNs are still a minor population in 
our sample, even at the faint end. The CTK 
AGNs identified above comprise 19 ± 6%, 9 ± 3%, 
3.6 ± 0.5%, and 1.9 ± 0.3% for the logSx ranges 
of (-16.1, -15.5), (-15.5, -15.0), (-15.0, -14.5), 
and (-14.5,-14.0) respectively, where the value 
of -16.1 corresponds to the faintest source in our 
sample. These fractions are consistent with the 
model prediction by U14 for /ctk = 1 (see Fig. 
17 of U14). 


3.11. Sample Summary 

In each sample except for Swift BAT, we have 
converted the 2-10 keV fluxes of all the sources 
using the energy conversion factor (ECF) assum¬ 
ing a power-law spectrum with a photon index of 
r = 1.7. To represent fluxes calculated in this 
way, we use the symbol Sx.obs, or S'x/,];/ when we 
want to emphasize that we use the ECF based on 
r = 1.7. We also define the observed luminosity 

Lx, obs = 47rd L (2) 2 Sx,obs, (1) 

where g?l(z) is the luminosity distance at the red- 
shift z. Implementation of the K-corrections and 
absorption corrections using realistic AGN spectra 
is discussed in Sect. 14.31 In our further analysis, 
we are interested in X-ray sources with intrinsic 
luminosities of log Lx > 42.0 and we use, in effect, 
X-ray point sources with log Lx,obs > 41.5. Thus 
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our analysis is strongly de-selecting non-AGN X- 
ray sources, the number density of which is lower 
than AGNs by a t least an order of magnitude at 
log L x > 4 1. 5 fPersic et all 12004 Ranal li et al 


2005 


Ptak et all 120071) . Our sample is plotted in 


the z — log Lx, obs space in Fig. [2j 

The properties of each sample used in our anal¬ 
ysis is summarized in Table [T] where the name 
of the survey, the survey solid angle, the limiting 
flux corresponding to 20% of the maximum solid 
angle, the number of extragalactic point sources 
with spectroscopic redshift, those with photomet¬ 
ric redshift, and identification completeness are 
listed for each survey. The numbers in this table 
do not include CTK AGNs that have been dis¬ 
cussed in the previous section nor those below our 
luminosity cut. 


3.12. Incompleteness Correction and Num¬ 
ber counts 


mulative number counts S^ 4ohs N(> Sx) curves 
in Fig. EJb), where Sxi4,obs is the 2-10 keV 
observed flux (see above) measured in units of 
10 _14 ergs _1 cm“ 2 , for each of our identified ex¬ 
tragalactic X-ray point source sample. The cumu¬ 
lative count is calculated by: 

N(> Sx, obs) = ^2 l/fi c (Sx,obs ,i), (2) 

-Sx ,obs,i >-Sx 

where f2 c (-Sx,obs,i) is the corrected survey area at 
the limiting flux of the i-th source -Sx, 0 bs,i- One 
sigma errors, calculated by 

<7 [N(> -Sx,obs)] 2 = 5] [lM(Sx,obs,0] 2 , 

Sx ,obs, >>-Sx ,obs 

(3) 

are shown. For visibility, we do not plot data 
points corresponding to the four brightest object 
of the survey. 


While the identification completeness of our se¬ 
lected samples are at least 90%, we further ac¬ 
count for unidentified sources as follows. For sam¬ 
ples with five or less unidentified sources, we have 
multiplied the survey area in the flux-limit — sur¬ 
vey area curve by the identification completeness 
of the sample over the whole flux range, assum¬ 
ing that the redshift (or the source type) proba¬ 
bility distribution of an unidentified source is the 
same as the redshift/source type distribution of 
the identified sources in the same sample. 

If there are more than five unidentified sources, 
we have made a similar adjustment of the survey 
area in flux blocks. We divide the X-ray sources 
into flux blocks from the brightest to faintest, in 
such a way that a new block starts after the fifth 
unidentified source. The last (faintest) block may 
contain one to five unidentified sources. For each 
block, we have adjusted the survey area by mul¬ 
tiplying by the completeness level of the block. 
The assumption is that the redshift/source type 
probability distribution is the same as that of the 
identified source at similar X-ray fluxes. 

The resulting incompleteness-corrected limiting- 
flux curve of each survey as well as that for the to¬ 
tal of all surveys are plotted in Fig. [3|a). In order 
to verify the integrity and consistency of our sam¬ 
ples and the level of systematic errors and/or the 
effects of cosmic variance, we plot the “tilted” cu¬ 


4. The XLF calculations 


4.1. Parametric Modeling with maximum 
likelihood fitting 


We follow the same procedu r e as our previous 
work ( Mivaii et al. 2000a . 2001 : Ileda et al. 2003, 


201411 and determine the best-fit parameters of pa¬ 


rameterized model of the XLF. 


As our maximum-likelihood estimator, we use 
either one of the following two forms: 


C = -2 ln Wi log Lxu Zi )] + 

2ff N (log L x , z)d log L x dz 


£ = -2 £> 


Af(log Lxi,Zj) _ 

J' J'N(log Lx,z)d log L x dz > 


( 4 ) 

( 5 ) 


where i goes through each AGN in the sample and 
N(\og Lx,z) (N-function) is the expected number 
density of AGNs in the sample per logarithmic lu¬ 
minosity per redshift, calculated from a parame¬ 
terized analytic model of the XLF: 


iV(log L x , z) = 

dL d Qg \ x d A {zf (1 + *) 2 &(*) • A(S X , obs), (6) 

where d A (z) is the angular diameter distance, % 
is the radial comoving distance (e.g., Sect 4.3 of 
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Table 1 

The X-ray samples 


Survey 

Max. Solid Angle 
[deg 2 ] 

lo gSx m obs a 

[cgs] 

AT b 
4V S Z 

iV pz c 

Completeness d 
% 

Swift BAT 

3.9 x 10 4 

-11.4 

274 

0 

98 

ASCA LSS+MSS 

91. 

-12.5 

94 

0 

97 

HBSS 

25. 

-12.6 

62 

0 

97 

Hellas2XMM 

3.8 

-13.5 

115 

0 

90 

XMS 

1.38 

-13.7 

128 

0 

91 

CLANS 

0.49 

-14.5 

183 

63 

90 

CLASXS 

0.32 

-14.4 

96 

23 

90 

XC-COSMOS 

2.17 

-14.6 

940 

424 

98 

LH-XMM 

0.20 

-14.7 

155 

102 

98 

CDF-N 

0.125 

-15.5 

118 

54 

90 

CDF-S 

0.128 

-15.8 

254 

85 

97 


a The limiting flux corresponding to 20% of the maximum solid angle of each survey. The 
fluxes are in the 2-10 keV band. 

b The number of extragalactic X-ray point sources with spectroscopic redshifts. 

c The number of extragalactic X-ray point sources with photometric redshifts, but with 
no spectroscopic redshifts. 

d Identification completeness. 




IQ-16 IQ-16 IQ-14 IQ-13 IQ-13 IQ-11 


Sx,o bs (2-10 keV) 


Sx.ob, (2-10 keV) 


Fig. 3.— (a) The completeness-corrected survey area curve as a function of the 2-10 keV flux limit. A 
photon index of T = 1.7 is assumed for all count-rate to flux conversions. The total area as well as those for 
individual surveys are shown as labeled, (b) The “tilted” cumulative number counts ob s N(> 5x,obs) for 
individual samples with lcr errors. The multiplication factor obs is applied to the cumulative number 
to make the Euclidean slope flat. The meaning of line colors and styles corresponds to those in the the flux 
limit - survey area curves in panel (a). For visibility, the data points corresponding to the 4 brightest objects 
in each sample are not displayed. (Color versions of these figures are available in the electronic version.) 
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Schneider 20061) . and A(Sx,obs) is the survey area 


as a function of limiting observed X-ray flux Sx,obs 
as shown in Fig. Ha). For a power-law spectrum 
with a photon index of T, 


ized with f p(z)dz = 1). The X-ray luminosity of 
the i — th object L Xj i(z) is now a function of z. 
The PDFs can be incorporated into Eqs. |4]&[5]by 
replacing A(log L xi ,Zi) by: 


<Sx,obs = Tx/47rd L (z) 2 (l + z) 2 r , (7) 


c?l(~) is the luminosity distance. 


The minimization of C with respect to model 
parameters gives the best-fit model. With Eq. 
H employed in our previous works, the parame¬ 
ter that represents the global normalization can¬ 
not be a fitting parameter, since within the esti¬ 
mator itself, the model number density is norrnal- 
ize d. On the othe r han d, wit h Eq. H employed 
by iMarshall et al.l (e.g., Fl983h . one can treat the 
normalization as a fitting parameter. 


We use the MINUIT package ( Jamesl[l994l) dis¬ 
tributed as a part of the CERN program library for 
the minimization procedure and parameter error 
search. The MINUIT command “MINOS” gives 
errors taking parameter correlations into account. 
In cases of global expressions (Sect. [5]), “MINOS” 
fails to give errors for many parameters if we use 
Eq. |4] In this case, we use Eq. [5] and the normal¬ 
ization and its errors are calculated independent of 
the fitting process as follows. The model normal¬ 
ization A (as a generic symbol of various normal¬ 
izations that appears in the following subsections) 
is determined, such that the total number of ex¬ 
pected objects is equal to the number of AGNs in 
the sample (A obs ). The estimated la confidence 
error for A is taken to be A/V A obs and does not 
include correlations of errors with other parame¬ 
ters. In our results, we indicate which method has 
been used for each fit. 


4.2. Photometric redshift full probability 
density distribution 

For the CLASXS, CLANS, COSMOS, LH, and 
Chandra Deep fields North and South datasets, 
significant fractions of redshifts rely on photomet¬ 
ric redshift determinations. Especially for fainter 
optical sources, which are systematically selected 
against for spectroscopic measurements, the error 
in the primary peak and catastrophic failure rate 
increase rapidly. Thus we include the probability 
density distribution (PDF) in redshift space for 
photometric redshifts, whenever available. 

Let the PDF of the i-th object be Pi{z) (normal- 


N (log L xi {z),z))pi(z)dz. 


( 8 ) 


We use the PDFs deriv e d by ISalvato et al.l 
( 2011 1 and Fotopoulou et al.l ( 2012l l for the COS¬ 
MOS and LH samples respectively. They give the 
PDFs in bins of A z = 0.01 in the 0 < z < 7 range. 
To reduce the computational time and at the same 
time not to sacrifice the accuracy, we rebinned the 
PDFs into Az/(1 + z) ~ 0.03 bins and then ne¬ 
glected the new bins that have a probability of 
less than 2%. For the CLASXS, CLANS, and 
CDF samples, no PDF information was available 
in the public domain at the ti me of ou r analy ¬ 


However, for the CDF-S, Xue et ahl ( 201ll) 


sis. 

show the first and second peaks of the photomet¬ 
ric redshift PDF derived by Luo et al. ( 2010h for 
some of their photometric redshifts. For those, we 
assumed that the PDF is the sum of two delta 
functions centered at these peaks, each of which 
has a weight of 0.5. The effects of considering the 
PDFs are discussed in Sect. 17.11 


4.3. Using Realistic AGN Spectra 

Our sample contains AGNs at various absorp¬ 
tion levels. Thus we have considered realistic ab¬ 
sorbed spectra in our XLF calculati ons, follow¬ 


ing the approach by U14 (see also iMivaii et al 


2000bl : Ileda et al. 2003 : Ebrero et al. 2009). In 
this work, we do not estimate the Ah on individ¬ 
ual AGNs nor derive the Ah function from our 
own dataset. Instead, we include in our model 
fitting the latest luminosity and redshift depen¬ 
dent Ah function, /(Lx, z; log Ah), from t he lit¬ 


eratu re and AGN template spectrum. Ueda et al 


( 2014| ) derived the refined Ah function and its evo¬ 
lution, where the evolution of absorbed AGN frac¬ 


tion is based on that of Has ingerj (2008). Writing 
down the full form of the Ah function is beyond 
the scope of this paper and the reader is directed 
to citetuedal4. The Ah function is normalized 
as J 2Q f(Lx, z\ log An)c?log Ah = 1 (unabsorbed 
AGNs with log Ah < 20.0 are included in the 
20 < log Ah < 21 bin). At z = 0, the Ah 
function of U14, which is based on the 9-month 
Swift BAT, AMSS and SXDS surveys, is in good 
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agreement with that of Burlon et ah (201l|), which 
is based on the three-year Swift BAT survey. 

Also U14 generated a set of sophisticated tem¬ 
plate X-ray spectra of AGNs, which include (i) 
an underlying power-law continuum with a high 
energy cutoff at E c = 300 keV, (ii) a reflection 
component, (iii) a scattering component from sur¬ 
rounding gas outside the torus, corresponding to 
1% with a torus opening solid angle is 27 t, and 
an intrinsic absorption with Compton-scattering 
derived from Monte-Carlo simulations. For our 
purpose, we take the spectral energy distribu¬ 
tion (SED) templates in Fig. 7. of U14 for 
log iV H [cm" 2 ] = 20.5,21.5,22.5, and 23.5, but ig¬ 
nored the luminosity and redshift dependence of 
the SED. In these template spectra, T = 1.94 and 
r = 1.84 are assumed for the primary power-law 
continuum for log Ah [cm -2 ] < 22 and > 22 AGNs 
respectively. 

In this section, we refine our definitions. The 
symbol Lx represents an intrinsic rest-frame 2-10 
keV luminosity of the primary power-law, which 
does not include the reflection component. This is 
because including the reflection component dou¬ 
ble counts the intrinsic isotropic luminosity of 
the AGN. The reflection component is typically 
~ 15% of the primary power-law in the 2-10 keV 
band. The symbol d<h/d\ogLx is the comoving 
number density of AGNs per dex in Lx for all 
Compton-thin (CTN) AGNs (log Ah < 24.0). 

To feed realistic spectra into the maximum like¬ 
lihood procedure, the A function, Eq. [G] should be 
a function of observable quantities for each sam¬ 
ple. We define the “observed” 2-10 keV luminosity 
defined in Eq. [T] 

In this case, clearly distinguishing between ob¬ 
served and intrinsic X-ray luminosities, Eq. [6] can 
be rewritten as: 


A (log Lx,obs,z)d log L x ,ohsdz = 
fw (dWT^ d l °gL x ) d A (z) 2 (1 + z) 2 ^{z)dz 
xf(L x , z; AH)A[5'x,obs(Ax, A H , z)]d\og A H , (9) 

where <Sx,obs is now a function of Ah in addi¬ 
tion to Lx and z. The integral is over log Ah 
and the intrinsic luminosity Lx is a unique 
function of Lx,obs, Ah, and z. The ratio of 
d log Lx,obs/d log Lx becomes non-unity only if 
the spectrum changes with Lx for a given Ah, 


which is the case for the Swift BAT sample (see 
below). 

In each template spectrum, with an intrinsic 
2-10 keV luminosity of the power-law continuum 
Lx, we can calculate the observed luminosity at 
2(1 + z)~ 10(1 + z) keV The 

true 2-10 keV flux from this object, under this 
template spectrum, is now: 

Sx,obs = Lx/*ndb(z) 2 • X ’ tm ^ % (10) 


The difference between S^bs ^xjobs 7 ( or the 
catalogued flux with an ECF calculated with T = 
1.7 or our assumed spectrum for the ECF) is often 
neglected in the literature. In reality, it becomes 
important for highly absorbed AGNs (logA H > 
23) at low redshifts \z < 1) in “2-10 keV” surveys, 
because the effective areas of the instruments used 
in generating our sample other than Swift BAT, 
drop rapidly as we go from 2 keV to 10 keV. Thus 
the ECF increases rapidly as the spectrum gets 
harder. 


An ideal approach would be to use separate re¬ 
sponse curve and countrate-survey area relation 
for each survey, as U14 did, since each survey 
use different instruments and/or different observa¬ 
tional epochs, for which the calibration may vary. 
However, a limitation of the software that we use 
for this work is that it can only accept one flux- 
limit - solid angle curve, as the software is a legacy 
from a series of single-ban d soft X -ray luminosity 
work Mivaii et all ( 2000al 2001 k Hasinger et al.l 
(12005 ). While this reduces the accuracy of the es¬ 
timate of 5*™ bs , the calculation speed gained by 
this simplification has allowed us to explore an¬ 
alytical forms for various sub-sets and cases (see 
below). We find that the relative ECFs between 
those for absorbed AGN spectra and for T = 1.7 
is not very sensitive to which instrument we use. 
For example, the ratios of the ECFs for an ex¬ 
tremely hard object (T = 1.7 with log Anh = 23.5 
at z=0.2) to that for an unabsorbed T = 1.7 spec¬ 
trum are 2.0, 2.1, and 1.7 for the XMM-Newton 
PN+2MOS (countrate at 2-8 kV), the Chandra, 
ACIS-I Cycle 8 (countrate at 2-7 keV) and the 
ASCA GIS (CR at 2-10 keV) respectively. Here we 
use the response for ACIS-I for the conversion be¬ 
tween 5‘( 0 u bs to Sr-'J. A caveat is that this causes 
some inaccuracy for the HBSS (CR is at 4.5-7.5 
keV) and XMM-Newton LH surveys (CR is at 2- 
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4.5 keV). However, the impact is relatively minor, 
since these latter surveys contain relatively small 
fraction of sources in the relevant flux range, i.e., 
the HBSS is overwhelmed by the ASCA MSS+LSS 
and the LH sample by the COSMOS sample in 
their respective flux coverages. 

For the Swift BAT sample, the selections are 
made using the 15-55 keV range, which is very lit¬ 
tle affected by the absorption up to log Ah = 24. 
Despite the limitation of our software described 
above, the Swift BAT sample can be treated sep¬ 
arately because the flux regime covered by this 
sample does not overlap with any other. We have 
made the conversion between the 2-10 keV and the 
15-55 keV bands as follows. 


In the comparison between Swift 14-194 keV 
and de-abs o rbed MAXI 2-10 keV luminosities, 
Ueda et al. ( 20111 ) pointed out that the effective 


power-law index between these bands varies with 
luminosity from Ti = 1.7 to T u = 2.0. We have 
confirmed this by comparing the 14 -55 keV lumi¬ 
nosities from the ( Aiello et al. 2012h catalog with 
4- 10 keV luminositie s in the latest MAXI catalog 
of lHiroi et al. (2013). To approximately represent 
the change, we use a smoothly varying luminosity- 
dependent effective spectral index 


T eff = Ti + 0.5 erf [(log Lx - 44.0)/2.0](T u - TO 

( 11 ) 

for the conversion between the 2 10 keV and 15-55 
keV fluxes, where erf is the error function. U14 hy¬ 
pothesize that the main cause of this slope change 
is the difference between type 1 and type 2 AGNs 
in the strength of component reflected by the torus 
(i.e., the torus tend to be thicker in type 2 AGNs) 
in combination with the relation that the type 2 
fraction decreases with luminosity. Because of the 
definition of Lx in this section, the 2-10 keV flux 
is calculated based on (1 + 0.15)Lx following the 
assumption that 15% of intrinsic 2-10 keV lumi¬ 
nosity is added by the reflection component. A 
small K-correction has been made based on T e ff. 
In order to use our flux limit - survey area curve, 
we have back-converted from the 15-55 keV flux 
to observed flux S£=b s 7 us i n § our reference slope 

r = i.7. 

As a technical note, we comment that U14 chose 
to use the count rate as an independent variable 
to the N-function instead of log Lx,obs- 


4.4. The A 0 b s /A m( ji method for the Binned 
XLF 


A less biased way of prese nting the binned 
XLF than the classical ]T)1/T4 (lAvini fe Bahcall 


198011 is the Anhs /Amdi estim ator descri bed in 
Mivaii et al. ( 200ll ) (see also lLa Franca fc Cristiani 


19971 ). The basic procedure is: 


1. Divide the combined sample into several red- 
shift shells. For each redshift shell, fit the 
AGN XLF with a smooth analytical func¬ 
tion obtained with the maximum-likelihood 
fit as described above. 

2. For each redshift shell, count the number 
of objects that fall into each luminosity bin 
to obtain the observed number of objects 
(Nobs)- 

3. For each luminosity bin, evaluate the ana¬ 
lytical fit at the central luminosity/redshift 
(d4> mdl /fflogL x ). 

4. Calculate the predicted number of AGNs in 
the bin (A md i). 

5. The final result is: 

d<b d4> mdl 

ttV = • N obs /A mdl . (12) 

d log L x d log L x 


We note that, in the case that the XLF model 
is constant within the bin, the A 0 b s /N md i esti¬ 
mator is mathematically equival ent to that pro¬ 
posed by [Page & Carrera (12000 ). In the case of 


taking absorption into account with the method 
described above in Sect. roi and binning is by 
log Lx defined there, deriving A 0 b s is not straight¬ 
forward, since we cannot determine log Lx indi¬ 
vidually from single-band fluxes. The Ah values 
derived from hardness ratios are subject to large 
errors and not necessarily available for all objects 
in the sample. On the other hand, it is straight¬ 
forward to calculate from the model probability 
distribution of Ah for given redshift and observed 
flux: p(z, S'x, obs ; log Ah), which can be normalized 
as: 


c24 


p(2,Sx,ob S ;logA H )dlog A h = 1. (13) 


/20 


This is different from the Ah function /(Lx, z; log Ah), 
which is the distribution of Ah for given redshift 


15 










































and intrinsic luminosity. For our purpose, for the 
*~th object with a given (zi, S x ,obs,i), we use the 
probability distribution of the intrinsic luminosity: 

Pi (log L x [(zi, S x , obM , N h )]) d log L x = 

p(zu S K>ohs ,p,logNji)dlogN H . (14) 


If we would like to calculate N 0 b s in the 
redshift-intrinsic luminosity bin (z\,Z 2 ', log L\, logl^): 


(■logL2 


N, 


obs 


= ^2 I Pi(}ogLx)dlogL x , 

JlogLi 


(15) 


where i is over objects with redshifts between z\ 
and Z 2 - In practice, we calculate -/V 0 b s as weighted 
sums over log Ah = 20.5,21.5,22.5 and 23.5 tem¬ 
plates. Similarly, we include the PDFs of photo¬ 
metric redshifts into N 0 b s estimations by adding 
the portion of the PDF of each objects that falls 
into the redshift bin, as has been made also by 
Aird et al.1 (|2010h . 

One limitation is that the errors assuming Pois¬ 
son statistics are not valid if A 0 b s is a sum of 
weights, unlike the cases of our previo us works 
( Mivaii et al. 20011 lHasinger et al.1 20051 U14), in 
which N ohs were always integers. However, we 
use the appro ximate lcr Poisson errors from the 
Oehrelsi (1986) Eqs. (7) and (12) (S = 1) for the 
upper and lower lcr errors respectively. Thus, the 
errors are only approximate. Also the errors in 
different bins are not independent of one another. 


5. Global XLF Expression 


A * is strongly coupled with the break luminosity 
Lx,*-, and therefore it is not possible to estimate 
its error accurately with our maximum-likelihood 
fitting procedure (see below). 

Following the most successful analytical form 
of the evolution in the literature, we use the 
luminosity-dependent density evolution: 


d $(L x ,z) d$(L x ,0) . r . 

d log Lx = d log L X 


(17) 


where e d (z,L x ) is the density evolution factor 
normalized to z=0. We describe it with a three- 
segment power-law model: 


e d (z,L x ) = 


(1 + z)p 1 

e d (*bi,£x)- (iTT^) 

/ \P3 

e d {zb 2 ,L x ) ■ 


(Z < Zbl) 

(Zbl < z < Zb2) 
(z > Z b2 ). 


(18) 

As in the previous works, the first break po¬ 
sitions z p i may depend on luminosity at the low 
luminosity regime: 


Zhi{L x ) 


~b,o(ix/ix,b)“ (L X < L x ,b) 

~b,o {L x > L Xj b) 


(19) 

We also include the luminosity dependence of 
the evolution indices pi and p 2 : 


Pi{L x ) = pi, 44 + Pi (log L x - 44) (20) 
P2{L X ) = P2,44 + 02 (log L x — 44). (21) 


5.1. The Luminosity-dependent Density 
Evolution Model 


For the analytical expression of the XLF at z=0, 
we use the smoothed two power-law formula. 


d$(L x ,0) 

d log L x 



The normalization A is the XLF value at 
(log Lx, z) = (44,0). There is also a convention 
to use A* (2nd line in Eq. fl6l) as a normaliza¬ 
tion parameter. We primarily use A|^°, because 


5.2. Global LDDE Results 

The best-fit parameters of our LDDE model 
and lcr errors, corresponding to AC = 1, are 
shown in Table [2] For the global expression, the 
fits using the likelihood function Eq. [5] give much 
more stable error search results than those using 
Eq. I As stated above, minimizing C in Eq. 
[5] cannot determine the normalization. Thus the 
A 44 0 value and its error are determined by scal¬ 
ing it such that it gives the total number of objects 
(Atot) for best-fit values of other parameters. The 
fractional error of l/V^tot is used for its error. 
Thus the errors of A 44 0 in Table [2] do not contain 
correlation of errors with other parameters. The 
traditional normalization AJ =0 (see Eq. m is also 
shown without an error for reference. 
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The parameters for the fit in which the real¬ 
istic AGN spectra, .ZVh function (Sect. 14.31) . and 
the PDFs of photometric redshifts (Sect. 14.21) are 
fully taken into account are listed under “Full” in 
Table [2] For our final results, we show the case 
where 2 bi ,44 is fixed to 1.1 and Zbi ,44 = 2.7, since 
fitting with all parameters free causes failures in 
the minimization processes. The values of these 
fixed parameters have been decided based on the 
results of the fits in luminosity-class divided sam¬ 
ples, which is explained later in Sect. 16.21 Since 
the fits in each luminosity class also use informa¬ 
tion on the global expression, thus the process has 
been iterative, i.e., we decide on fixed parameters 
of the global expression referring to the results of 
the luminosity-class divided sample based on the 
global expression in the previous step. 

6. XLFs in redshift shells and luminosity 
classes 

6 . 0 . 1 . Adaptive binning 

In this section, we present the binned XLF, as 
well as analytical formula of the XLF divided into 
redshift shells. We also present the number den¬ 
sity curve divided in luminosity classes. The bin¬ 
ning has been decided on a adaptive-basis, so that 
the sizes of error bars and bin-sizes are optimal for 
the plots to be visible and informative. 

First, we divide the sample into 11 redshift 
shells in the same way as those used by U14, ex¬ 
cept that the two highest redshift shells by U14 
(z > 3.0) are combined into a single one for obtain¬ 
ing a sufficient number of objects for a separate 
two power-law fit. For each redshift shell, we con¬ 
struct luminosity bins such that each contains 20 
objects by default, with minimum and maximum 
bin sizes of A log Lx = 0.125 and 0.5 respectively. 
In case the most luminous bin contains less than 
20 objects, the border between the last two bins is 
adjusted to contain approximately equal numbers 
of objects within the bin size constraints. 

For the luminosity division, the range 42.0 < 
log Lx < 46.0 is divided into four A log Lx = 1.0 
classes. Each luminosity class is divided into red¬ 
shift bins such that each bin has 20 objects (using 
the best fit photometric redshifts for those with¬ 
out spectroscopic redshift) with forced minimum 
and maximum sizes of Aln(l + z) = 0.07 and 0.5 
respectively. If the highest redshift bin has less 


than 20 objects, the border between the highest 
two redshift bins are adjusted in the same way as 
the case of luminosity binning. Due to the small 
number of objects, the default number of objects 
per redshift bin for the 45.0 < log Lx < 46.0 class 
has been set to 7, with the same minimum and 
maximum bin size constraints. 


6.1. Analytical form for each redshift shell 

For the analytical expression of the XLF in each 
redshift shell, we use the smoothed two power law 
formula. Because the redshift shells have finite 
widths, the fit results depend on the evolution of 
the SXLF within them: 

d ( f> (L x , z) 
d log L x 

x ea (z,L x ), (22) 

where e A (z,Lx) is the density evolution factor. 
While the final results are insensitive to the de¬ 
tailed behavior of e<i(z, Lx) within the shell at 
most locations in the (Lx, z) space, we have 
taken our best-estimate by using the luminosity- 
dependent density evolution (LDDE) model de¬ 
rived above in Sect. [5] The luminosity range of 
the fit is from log Lx,obs=41.5 to 46.0. 

The best fit parameters of Eq. [55] for each red¬ 
shift shell are shown in Tabled The normalization 
is defined by: 

A z= Zc _ d $ (L x = 10 44 erg s~ 1 ,z = z c ) 

44 d log L x ’ ( } 

where z c is the central redshift of the shell, which is 
defined by ~ c + 1 = i/(l + z m in)(l + Anax), where 
-min and Zmax are the minimum and maximum 
redshifts of the bin. One sigma parameter errors 
are calculated in the same way as in Sect. 15.21 
In this section, we have been able to obtain la 
errors of most parameters with Eq. [U and unlike 
the case of Sect. 15.21 the errors of A aa Zc contain 
the correlations among parameter errors. We also 
show the traditional normalization A z=Zc without 
errors. 

Table Q] tabulates the binned Aobs/Andi- 
estimated XLF for the full treatment case, along 
with the observed number of AGNs (N 0 b s ) and 
the final estimated values XLF value at the cen¬ 
ter of each bin. The full XLFs in the 11 redshift 
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Table 2 

Best-fit LDDE parameters for Global Expressions 


Parameter* 1 Full No pz-PDZ F = 1.8 


z=0 XLF Parameters 13 


AZ = Oc 

^44 

(9.3 ±0.2) x 10“ 7 

(8.7 ±0.2) x 10" 7 

(9.7 ±0.2) X 10- 7 

A z= - Od 

1.56 x 10 _ ® 

1.34 x lO” 6 

1.60 x 10“® 

log ±x,» 

44.04 ± 0.08 

44.06 ± 0.06 

44.05 ± 0.06 

71 

1.17 ±0.05 

1.17 ±0.05 

1.09 ±0.04 

72 

<o O 

o d 
+ 1 
o 

GO 

CN 

2 90 +0 ' 08 
z.yu_o 12 

2.76 ±0.08 

Evolution Parameters 13 

Pi,44 

5.29 ±0.11 

5.35 ±0.11 

4.86 ±0.11 

Zbl,44 

1.1* 

1.1* 

1.1* 

P2,44 

-0.35 ±0.14 

0.02 ± 0.13 

0.02 ±0.14 

z b2 

2.7* 

2.7* 

2.7* 

p3 44 

-5.6* 

-5.6* 

-5.6* 

a 

0.18 ± .03 

0-16t ,oi 

0 11 + ' 01 
u - i± -,02 

logLx.b 

44.5* 

44.5* 

44.5* 

Pi 

i O+0.3 

1 - z -0.2 

i 9+0.2 

1 -0.3 

1.6 ±0.2 

02 

i c+0.5 

1 -°-0.3 

2 (V°' 2 

4- u _o,4 

1.7 ± 0.2 


a Parameters that have been fixed during the fit are labeled by ’(*)’. Units 
- A (with any sub and superscripts): h. 3 0 Mpc -3 dex — L (with any sub and 
superscripts): h^ergs -1 . 

b The 68% confidence range for one parameter (A C < 1) with correlations 
among parameter errors, except for A^°. 

c The lcr errors of AJJ 0 is for all other parameters fixed. 

d The traditional normalization defined as the twice of the XLF value at the 
break luminosity at z = 0 (L x ,*). 


Table 3 

Smoothed 2PL parameters for each redshift shell 


^-range 

2c 

A z — Zc 
2^44 

a 

b 

Af 

-Zc a b 

log L« — Zc b 

71 a 

72 a 

0.015 - 

0.200 

0.104 

i.zu_ 0 13 

X 

10 -6 

8.97 

X 

10-® 

43.70++ 

0.90+;!! 

9 rq+0.23 
z.ao-o.19 

0.200 - 

0.400 

0.296 

4 7Q+ 0 - 65 
' y —0.61 

X 

10 _ ® 

3.91 

X 

10-® 

44.22«; 2 i 

L05+;°® 

2.90!°;| 7 

0.400 - 

0.600 

0.497 

' 1.38 

X 

10~® 

1.67 

X 

10- 5 

43.98++ 

0.83+;!! 

2.70l°;« 

0.600 - 

0.800 

0.697 

i 14+O.17 
1 * 14 -0.14 

X 

10- 5 

2.36 

X 

10" 4 

43.31++ 

0 41 + 0 - 24 

u - 41 _0.30 

1-86++ 

0.800 - 

1.000 

0.897 

q 1 O+0.37 
°- ±z —0.35 

X 

10~ 5 

8.91 

X 

10- 5 

43.90++ 

0 4Q+ 013 
u -^ u -0.15 

9 ci +0.27 
z>o± -0.22 

1.000 - 

1.200 

1.098 

0 4-7+0.35 
_q.33 

X 

10 _s 

1.06 

X 

10- 4 

43.85+;!! 

n iq+ 016 
u - iy -0.19 

2.06!°;!® 

1.200 - 

1.600 

1.392 

4.42+31 

X 

10~ 5 

3.04 

X 

10- 5 

44.43+;°® 

0-47+;° 7 

9 09+O.32 

Z - OZ -0.27 

1.600 - 

2.000 

1.793 

3-68++ 

X 

10- 5 

2.07 

X 

10" 5 

44.57+;!° 

0.48+;°° 

9 -70+0.38 

'°-0.32 

2.000 - 

2.400 

2.194 

4 9-7+0.50 
-0.43 

X 

10~ 5 

3.28 

X 

10- 5 

44.50+;!! 

0.27+;!! 

2.84!°;!! 

2.400 - 

3.000 

2.688 

q 90 + 0.48 
O.ZO-033 

X 

10" s 

1.71 

X 

10- 5 

44.63++ 

0.46+;!! 

3-1+aeo 

3.000 - 

5.800 

4.215 

4 q-7+0.90 
4,0 *-0.62 

X 

10~® 

1.08 

X 

10-® 

44 Q4"*"2*H9 
^• y ^_0.28 

0.65+i° 

4-38!+ 9 


a Errors show the 68% confidence range for one parameter (A C, < 1). 

b Units - A (with any sub/upper scripts): [h 3 0 Mpc -3 dex — *], L (with any sub/upper scripts): 

[ Vo 2er s S_1 1- 
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shells are plotted in Fig.[4]in separate panels. The 
best-fit smoothed two power-law model (2PL) and 
the best-fit LDDE model evaluated at z c are over¬ 
plotted to assess the goodness of the models. The 
best-fit two power law models at z c = 0.1 (low¬ 
est redshift bin) and at z c = 1.8, where the A 44 
is near the peak, are also shown in all panels for 
reference. 
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42 43 44 45 42 43 44 45 42 43 44 45 


log L x [h^f erg s" 1 ] 

Fig. 4.— The N a b s /N m di estimated intrinsic 2-10 keV X-ray luminosity function is plotted in each redshift 
shell as labeled. The 68% approximate Poisson error bars are shown in each data point. The solid line in 
each panel shows the redshift divided smoothed 2PL model at the central redshift (z c ). The best-fit LDDE 
model is shown with dotted lines. The best-fit smoothed two power law models for the lowest redshift shell 
(z c = 0.093) and at z c = 2.214 are shown with dashed lines in all other panels for reference. (Color version 
of this figure is available in the electronic version.) 
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Table 4 

Full binned XLF divided in redshift shells 


z-range 

2c 

log L x - 

range a 

log L x , c a 

-^obj 


d 

d$ a 

log Lx 

b 


0.015 - 

0.200 

0.104 

42.13 - 

42.48 

42.31 

11.4 

(1.24 

± 

0.24) 

X 

10“ 

-4 

0.015 - 

0.200 

0.104 

42.48 - 

42.61 

42.54 

5.2 

(5.85 

± 

1.76) 

X 

10“ 

-5 

0.015 - 

0.200 

0.104 

42.61 - 

42.73 

42.67 

15.9 

(1.02 

± 

0.23) 

X 

10" 

-4 

0.015 - 

0.200 

0.104 

42.73 - 

42.86 

42.79 

16.1 

(5.58 

± 

1.28) 

X 

10- 

-5 

0.015 - 

0.200 

0.104 

42.86 - 

42.98 

42.92 

20.0 

(4.68 

± 

0.89) 

X 

10“ 

-5 

0.015 - 

0.200 

0.104 

42.98 - 

43.11 

43.04 

26.1 

(3.61 

± 

0.65) 

X 

10- 

-5 

0.015 - 

0.200 

0.104 

43.11 - 

43.23 

43.17 

26.4 

(2.54 

± 

0.46) 

X 

10- 

-5 

0.015 - 

0.200 

0.104 

43.23 - 

43.36 

43.29 

26.6 

(1.54 

± 

0.26) 

X 

10“ 

-5 

0.015 - 

0.200 

0.104 

43.36 - 

43.48 

43.42 

28.7 

(1.20 

± 

0.20) 

X 

10" 

-5 

0.015 - 

0.200 

0.104 

43.48 - 

43.61 

43.54 

29.6 

(7.81 

± 

1.33) 

X 

10“ 

-6 

0.015 - 

0.200 

0.104 

43.61 - 

43.73 

43.67 

28.2 

(5.55 

± 

1.00) 

X 

10“ 

-6 

0.015 - 

0.200 

0.104 

43.73 - 

43.86 

43.79 

14.8 

(1.85 

± 

0.46) 

X 

io- 

-6 

0.015 - 

0.200 

0.104 

43.86 - 

43.98 

43.92 

24.4 

(2.25 

± 

0.45) 

X 

10“ 

-6 

0.015 - 

0.200 

0.104 

43.98 - 

44.11 

44.04 

18.9 

(1.08 

± 

0.25) 

X 

10“ 

-6 

0.015 - 

0.200 

0.104 

44.11 - 

44.23 

44.17 

12.5 

(5.18 

± 

1.45) 

X 

10“ 

-7 

0.015 - 

0.200 

0.104 

44.23 - 

44.42 

44.32 

9.2 

(1.44 

± 

0.48) 

X 

10“ 

-7 

0.015 - 

0.200 

0.104 

44.42 - 

44.55 

44.49 

10.0 

(1.34 

± 

0.43) 

X 

io- 

-7 

0.015 - 

0.200 

0.104 

44.55 - 

44.70 

44.63 

11.0 

(7.88 

± 

2.36) 

X 

10" 

-8 

0.015 - 

0.200 

0.104 

44.70 - 

45.10 

44.90 

7.0 

(7.51 

± 

2.85) 

X 

10“ 

-9 

0.015 - 

0.200 

0.104 

45.10 - 

45.35 

45.22 

2.0 

(2.53 

± 

1.77) 

X 

10“ 

-9 

0.200 - 

0.400 

0.296 

41.89 - 

42.07 

41.98 

5.9 

(2.48 

± 

0.62) 

X 

10“ 

-4 

0.200 - 

0.400 

0.296 

42.07 - 

42.30 

42.18 

18.4 

(3.50 

± 

0.53) 

X 

10“ 

-4 

0.200 - 

0.400 

0.296 

42.30 - 

42.42 

42.36 

13.3 

(3.11 

± 

0.53) 

X 

io- 

-4 

0.200 - 

0.400 

0.296 

42.42 - 

42.55 

42.49 

12.7 

(2.03 

± 

0.35) 

X 

10“ 

-4 

0.200 - 

0.400 

0.296 

42.55 - 

42.67 

42.61 

16.8 

(2.21 

± 

0.33) 

X 

io- 

-4 

0.200 - 

0.400 

0.296 

42.67 - 

42.80 

42.74 

14.3 

(1.36 

± 

0.20) 

X 

10“ 

-4 

0.200 - 

0.400 

0.296 

42.80 - 

42.92 

42.86 

12.2 

(9.95 

± 

1.59) 

X 

10“ 

-5 

0.200 - 

0.400 

0.296 

42.92 - 

43.05 

42.99 

10.1 

(6.07 

± 

1.03) 

X 

io- 

-5 

0.200 - 

0.400 

0.296 

43.05 - 

43.17 

43.11 

15.0 

(7.75 

± 

1.08) 

X 

10" 

-5 

0.200 - 

0.400 

0.296 

43.17 - 

43.30 

43.24 

11.7 

(4.67 

± 

0.75) 

X 

10- 

-5 

0.200 - 

0.400 

0.296 

43.30 - 

43.48 

43.39 

13.5 

(3.17 

± 

0.48) 

X 

10- 

-5 

0.200 - 

0.400 

0.296 

43.48 - 

43.60 

43.54 

4.9 

(1.25 

± 

0.33) 

X 

10- 

-5 

0.200 - 

0.400 

0.296 

43.60 - 

43.80 

43.70 

12.1 

(1.08 

± 

0.18) 

X 

10- 

-5 

0.200 - 

0.400 

0.296 

43.80 - 

43.92 

43.86 

8.4 

(6.79 

± 

1.56) 

X 

10“ 

-6 

0.200 - 

0.400 

0.296 

43.92 - 

44.05 

43.99 

11.4 

(5.31 

± 

1.12) 

X 

10“ 

-6 

0.200 - 

0.400 

0.296 

44.05 - 

44.23 

44.14 

11.2 

(2.49 

± 

0.52) 

X 

10- 

-6 

0.200 - 

0.400 

0.296 

44.23 - 

44.51 

44.37 

9.8 

(1.01 

± 

0.24) 

X 

10" 

-6 

0.200 - 

0.400 

0.296 

44.51 - 

45.01 

44.76 

4.3 

(1.28 

± 

0.51) 

X 

10- 

-7 

0.200 - 

0.400 

0.296 

45.01 - 

45.16 

45.08 

1.1 

(9.41 

± 

8.85) 

X 

10- 

-9 

0.200 - 

0.400 

0.296 

45.16 - 

45.56 

45.36 

2.0 

(1.28 

± 

0.90) 

X 

10- 

-9 

0.200 - 

0.400 

0.296 

45.56 - 

45.81 

45.69 

2.0 

(5.77 

± 

4.04) 

X 

10- 

10 

0.400 - 

0.600 

0.497 

41.90 - 

42.19 

42.04 

5.9 

(2.92 

± 

0.70) 

X 

10- 

-4 

0.400 - 

0.600 

0.497 

42.19 - 

42.31 

42.25 

6.3 

(3.58 

± 

0.86) 

X 

10- 

-4 

0.400 - 

0.600 

0.497 

42.31 - 

42.46 

42.39 

10.5 

(2.96 

± 

0.56) 

X 

10“ 

-4 

0.400 - 

0.600 

0.497 

42.46 - 

42.59 

42.53 

9.8 

(1.96 

± 

0.39) 

X 

10“ 

-4 

0.400 - 

0.600 

0.497 

42.59 - 

42.71 

42.65 

15.3 

(2.29 

± 

0.34) 

X 

10“ 

-4 

0.400 - 

0.600 

0.497 

42.71 - 

42.85 

42.78 

22.1 

(1.95 

± 

0.23) 

X 

10- 

-4 

0.400 - 

0.600 

0.497 

42.85 - 

42.98 

42.92 

22.7 

(1.53 

± 

0.18) 

X 

10- 

-4 
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Table 4 —Continued 


2-range 

2 c 

log L x - 

range a 

logL x ,o a 

-N°bj 


d 

d < E > 

log Ly 

_ a 

,b 

0.400 - 

0.600 

0.497 

42.98 - 

43.10 

43.04 

17.0 

( 9.81 

± 

1 . 28 ) 

X 

10 “ 

0.400 - 

0.600 

0.497 

43.10 - 

43.23 

43.17 

17.3 

( 7.55 

± 

0 . 98 ) 

X 

10 ” 

0.400 - 

0.600 

0.497 

43.23 - 

43.35 

43.29 

18.0 

( 6.69 

± 

0 . 87 ) 

X 

10 " 

0.400 - 

0.600 

0.497 

43.35 - 

43.48 

43.42 

16.7 

( 4.45 

± 

0 . 58 ) 

X 

10 - 

0.400 - 

0.600 

0.497 

43.48 - 

43.60 

43.54 

12.7 

( 2.85 

± 

0 . 43 ) 

X 

10 - 

0.400 - 

0.600 

0.497 

43.60 - 

43.73 

43.67 

11.8 

( 2.08 

± 

0 . 33 ) 

X 

io - 

0.400 - 

0.600 

0.497 

43.73 - 

43.88 

43.81 

11.8 

( 1.60 

± 

0 . 27 ) 

X 

io - 

0.400 - 

0.600 

0.497 

43.88 - 

44.25 

44.07 

13.2 

( 4.79 

± 

0 . 77 ) 

X 

io - ( 

0.400 - 

0.600 

0.497 

44.25 - 

44.41 

44.33 

7.9 

( 2.02 

± 

0 . 55 ) 

X 

io - 1 

0.400 - 

0.600 

0.497 

44.41 - 

44.53 

44.47 

6.5 

( 1.25 

± 

0 . 39 ) 

X 

io - 1 

0.400 - 

0.600 

0.497 

44.53 - 

44.77 

44.65 

2.3 

( 1.42 

± 

0 . 54 ) 

X 

io —' 

0.400 - 

0.600 

0.497 

44.77 - 

45.13 

44.95 

2.1 

( 6.14 

± 

3 . 50 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

41.88 - 

42.36 

42.12 

10.9 

( 4.06 

± 

0 . 69 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

42.36 - 

42.48 

42.42 

6.3 

( 3.83 

± 

0 . 84 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

42.48 - 

42.61 

42.54 

11.8 

( 4.05 

± 

0 . 69 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

42.61 - 

42.73 

42.67 

18.4 

( 4.20 

± 

0 . 55 ) 

X 

10 “ 

0.600 - 

0.800 

0.697 

42.73 - 

42.86 

42.79 

20.9 

( 2.99 

± 

0 . 39 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

42.86 - 

42.98 

42.92 

29.1 

( 3.26 

± 

0 . 33 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

42.98 - 

43.11 

43.04 

29.0 

( 1.96 

± 

0 . 20 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.11 - 

43.23 

43.17 

36.7 

( 1.83 

± 

0 . 17 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.23 - 

43.36 

43.29 

30.1 

( 1.08 

± 

0 . 11 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.36 - 

43.48 

43.42 

30.5 

( 9.83 

± 

0 . 95 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.48 - 

43.61 

43.54 

21.1 

( 5.09 

± 

0 . 56 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.61 - 

43.73 

43.67 

23.0 

( 4.68 

± 

0 . 51 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.73 - 

43.86 

43.79 

14.3 

( 2.08 

± 

0 . 29 ) 

X 

10 - 

0.600 - 

0.800 

0.697 

43.86 - 

43.98 

43.92 

13.3 

( 1.68 

± 

0 . 25 ) 

X 

10 “ 

0.600 - 

0.800 

0.697 

43.98 - 

44.13 

44.06 

9.5 

( 8.32 

± 

1 . 58 ) 

X 

10 -' 

0.600 - 

0.800 

0.697 

44.13 - 

44.52 

44.33 

12.5 

( 3.20 

± 

0 . 54 ) 

X 

10 -' 

0.600 - 

0.800 

0.697 

44.52 - 

44.65 

44.58 

5.0 

( 1.52 

± 

0 . 53 ) 

X 

10 -' 

0.600 - 

0.800 

0.697 

44.65 - 

45.08 

44.87 

15.5 

( 4.40 

± 

0 . 88 ) 

X 

10 —' 

0.600 - 

0.800 

0.697 

45.08 - 

45.33 

45.21 

0.2 

( 7.45 

± 

2 . 46 ) 

X 

10 “' 

0.800 - 

1.000 

0.897 

41.98 - 

42.48 

42.23 

3.9 

( 1.49 

± 

0 . 43 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

42.48 - 

42.73 

42.60 

12.8 

( 3.40 

± 

0 . 51 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

42.73 - 

42.85 

42.79 

11.5 

( 3.13 

± 

0 . 47 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

42.85 - 

42.98 

42.92 

15.5 

( 2.27 

± 

0 . 27 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

42.98 - 

43.10 

43.04 

22.2 

( 2.36 

± 

0 . 26 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.10 - 

43.23 

43.17 

25.3 

( 1.75 

± 

0 . 19 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.23 - 

43.35 

43.29 

25.4 

( 1.26 

± 

0 . 13 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.35 - 

43.48 

43.42 

41.2 

( 1.39 

± 

0 . 11 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.48 - 

43.60 

43.54 

39.3 

( 1.20 

± 

0 . 10 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.60 - 

43.73 

43.67 

39.0 

( 8.68 

± 

0 . 74 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.73 - 

43.85 

43.79 

28.5 

( 5.37 

± 

0 . 54 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.85 - 

43.98 

43.92 

28.2 

( 3.97 

± 

0 . 40 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

43.98 - 

44.10 

44.04 

16.5 

( 1.93 

± 

0 . 25 ) 

X 

10 - 

0.800 - 

1.000 

0.897 

44.10 - 

44.23 

44.17 

19.7 

( 1.69 

± 

0 . 22 ) 

X 

10 “ 

0.800 - 

1.000 

0.897 

44.23 - 

44.48 

44.35 

18.3 

( 6.69 

± 

0 . 94 ) 

X 

10 "' 

0.800 - 

1.000 

0.897 

44.48 - 

44.88 

44.68 

10.5 

( 1.50 

± 

0 . 33 ) 

X 

10 "' 

0.800 - 

1.000 

0.897 

44.88 - 

45.13 

45.01 

1.3 

( 6.39 

± 

3 . 90 ) 

X 

10 -' 
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Table 4 —Continued 


z -range 

2 c 

log L x - 

range a 

logL x , o a 

^obj 


d 

d < E > 

log Ly 

_ a 

,b 

1.000 - 

1.200 

1.098 

41.86 - 

42.36 

42.11 

1.8 

( 1.14 

± 

0 . 42 ) 

X 

10 - 

1.000 - 

1.200 

1.098 

42.36 - 

42.85 

42.60 

11.1 

( 2.15 

± 

0 . 37 ) 

X 

10 “ 

1.000 - 

1.200 

1.098 

42.85 - 

42.98 

42.92 

6.4 

( 1.76 

± 

0 . 35 ) 

X 

10 “ 

1.000 - 

1.200 

1.098 

42.98 - 

43.10 

43.04 

10.4 

( 1.80 

± 

0 . 29 ) 

X 

io - 

1.000 - 

1.200 

1.098 

43.10 - 

43.23 

43.17 

14.2 

( 1.37 

± 

0 . 18 ) 

X 

io - 

1.000 - 

1.200 

1.098 

43.23 - 

43.35 

43.29 

14.9 

( 1.12 

± 

0 . 15 ) 

X 

io - 

1.000 - 

1.200 

1.098 

43.35 - 

43.48 

43.42 

22.7 

( 1.04 

± 

0 . 11 ) 

X 

IO " 

1.000 - 

1.200 

1.098 

43.48 - 

43.60 

43.54 

30.4 

( 1.02 

± 

0 . 09 ) 

X 

IO " 

1.000 - 

1.200 

1.098 

43.60 - 

43.73 

43.67 

31.3 

( 7.51 

± 

0 . 68 ) 

X 

io - 

1.000 - 

1.200 

1.098 

43.73 - 

43.85 

43.79 

30.3 

( 6.65 

± 

0 . 63 ) 

X 

io - 

1.000 - 

1.200 

1.098 

43.85 - 

43.98 

43.92 

27.6 

( 4.53 

± 

0 . 45 ) 

X 

IO " 

1.000 - 

1.200 

1.098 

43.98 - 

44.10 

44.04 

24.6 

( 3.66 

± 

0 . 40 ) 

X 

IO " 

1.000 - 

1.200 

1.098 

44.10 - 

44.23 

44.17 

16.7 

( 1.77 

± 

0 . 23 ) 

X 

IO " 

1.000 - 

1.200 

1.098 

44.23 - 

44.35 

44.29 

12.6 

( 1.07 

± 

0 . 17 ) 

X 

io - 

1.000 - 

1.200 

1.098 

44.35 - 

44.57 

44.46 

11.8 

( 4.37 

± 

0 . 74 ) 

X 

io -' 

1.000 - 

1.200 

1.098 

44.57 - 

44.98 

44.78 

11.4 

( 1.69 

± 

0 . 34 ) 

X 

io -' 

1.000 - 

1.200 

1.098 

44.98 - 

45.34 

45.16 

4.3 

( 1.53 

± 

0 . 58 ) 

X 

10 —' 

1.000 - 

1.200 

1.098 

45.34 - 

45.59 

45.47 

4.0 

( 8.73 

± 

3 . 75 ) 

X 

10- 1 

1.200 - 

1.600 

1.392 

42.10 - 

42.60 

42.35 

7.0 

( 1.82 

± 

0 . 40 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

42.60 - 

42.86 

42.73 

6.3 

( 1.34 

± 

0 . 31 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

42.86 - 

43.04 

42.95 

6.7 

( 1.10 

± 

0 . 24 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.04 - 

43.16 

43.10 

10.6 

( 1.53 

± 

0 . 26 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.16 - 

43.29 

43.22 

14.9 

( 1.15 

± 

0 . 15 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.29 - 

43.41 

43.35 

23.1 

( 1.17 

± 

0 . 13 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.41 - 

43.54 

43.47 

26.9 

( 8.54 

± 

0 . 74 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.54 - 

43.66 

43.60 

37.6 

( 8.73 

± 

0 . 71 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.66 - 

43.79 

43.72 

49.2 

( 7.16 

± 

0 . 49 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.79 - 

43.91 

43.85 

50.6 

( 5.98 

± 

0 . 41 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

43.91 - 

44.04 

43.97 

58.3 

( 5.17 

± 

0 . 34 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

44.04 - 

44.16 

44.10 

44.0 

( 3.40 

± 

0 . 26 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

44.16 - 

44.29 

44.22 

39.1 

( 2.31 

± 

0 . 19 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

44.29 - 

44.41 

44.35 

38.2 

( 2.04 

± 

0 . 18 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

44.41 - 

44.54 

44.47 

25.8 

( 1.03 

± 

0 . 11 ) 

X 

10 - 

1.200 - 

1.600 

1.392 

44.54 - 

44.66 

44.60 

23.7 

( 8.19 

± 

1 . 06 ) 

X 

10 -' 

1.200 - 

1.600 

1.392 

44.66 - 

44.98 

44.82 

21.9 

( 2.17 

± 

0 . 26 ) 

X 

10 -' 

1.200 - 

1.600 

1.392 

44.98 - 

45.10 

45.04 

2.0 

( 4.94 

± 

2 . 03 ) 

X 

10 —' 

1.200 - 

1.600 

1.392 

45.10 - 

45.29 

45.19 

4.2 

( 3.90 

± 

1 . 44 ) 

X 

10 —' 

1.200 - 

1.600 

1.392 

45.29 - 

45.78 

45.53 

2.2 

( 1.48 

± 

0 . 77 ) 

X 

10 _! 

1.200 - 

1.600 

1.392 

45.78 - 

46.03 

45.90 

1.0 

( 7.86 

± 

7 . 00 ) 

X 

10 " ! 

1.600 - 

2.000 

1.793 

41.91 - 

42.41 

42.16 

1.3 

( 1.45 

± 

0 . 62 ) 

X 

10 - 

1.600 - 

2.000 

1.793 

42.41 - 

42.91 

42.66 

7.3 

( 1.38 

± 

0 . 30 ) 

X 

10 - 

1.600 - 

2.000 

1.793 

42.91 - 

43.21 

43.06 

8.0 

( 1.05 

± 

0 . 22 ) 

X 

10 - 

1.600 - 

2.000 

1.793 

43.21 - 

43.38 

43.29 

6.5 

( 7.14 

± 

1 . 71 ) 

X 

10 “ 

1.600 - 

2.000 

1.793 

43.38 - 

43.51 

43.44 

9.2 

( 6.95 

± 

1 . 11 ) 

X 

10 - 

1.600 - 

2.000 

1.793 

43.51 - 

43.63 

43.57 

14.8 

( 7.12 

± 

0 . 85 ) 

X 

10 - 

1.600 - 

2.000 

1.793 

43.63 - 

43.76 

43.69 

21.8 

( 6.48 

± 

0 . 63 ) 

X 

10 - 

1.600 - 

2.000 

1.793 

43.76 - 

43.88 

43.82 

22.2 

( 4.85 

± 

0 . 48 ) 

X 

10 “ 

1.600 - 

2.000 

1.793 

43.88 - 

44.01 

43.94 

38.1 

( 5.16 

± 

0 . 37 ) 

X 

10 “ 

1.600 - 

2.000 

1.793 

44.01 - 

44.13 

44.07 

30.5 

( 3.30 

± 

0 . 28 ) 

X 

10 - 
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Table 4 —Continued 


z -range 

2 c 

log Lx-range a 

logLx,c a 

-N°bj 



d<E> 

_ a 

,b 

d log Ly 


1.600 - 2.000 

1.793 

44.13 - 44.26 

44.19 

38.0 

(3.07 

± 

0.24) 

X 

10” 5 

1.600 - 2.000 

1.793 

44.26 - 44.38 

44.32 

32.7 

(2.30 

± 

0 .20) 

X 

10” 5 

1.600 - 2.000 

1.793 

44.38 - 44.51 

44.44 

25.3 

(1.35 

± 

0.14) 

X 

10” 5 

1.600 - 2.000 

1.793 

44.51 - 44.63 

44.57 

18.4 

(9.15 

± 

1 .10) 

X 

10 ”® 

1.600 - 2.000 

1.793 

44.63 - 44.80 

44.72 

16.2 

(4.14 

± 

0.54) 

X 

10 ”® 

1.600 - 2.000 

1.793 

44.80 - 44.93 

44.86 

8.3 

(2.14 

± 

0.45) 

X 

10 ”® 

1.600 - 2.000 

1.793 

44.93 - 45.05 

44.99 

3.2 

(8.10 

± 

2.59) 

X 

10- 7 

1.600 - 2.000 

1.793 

45.05 - 45.18 

45.11 

3.4 

(7.08 

± 

2.55) 

X 

10- 7 

1.600 - 2.000 

1.793 

45.18 - 45.61 

45.40 

2.1 

(7.76 

± 

2.87) 

X 

10” 8 

1.600 - 2.000 

1.793 

45.61 - 45.86 

45.74 

0.9 

(1.14 

± 

1 .01) 

X 

10” 8 

2.000 - 2.400 

2.194 

42.58 - 43.08 

42.83 

4.4 

(8.19 

± 

2.29) 

X 

10 ” s 

2.000 - 2.400 

2.194 

43.08 - 43.45 

43.26 

7.1 

(6.74 

± 

1.55) 

X 

10” 5 

2.000 - 2.400 

2.194 

43.45 - 43.63 

43.54 

6.9 

(5.42 

± 

1.08) 

X 

10 ” s 

2.000 - 2.400 

2.194 

43.63 - 43.76 

43.69 

9.7 

(5.26 

± 

0.74) 

X 

10” 5 

2.000 - 2.400 

2.194 

43.76 - 43.88 

43.82 

14.8 

(5.46 

± 

0 .66) 

X 

10” 5 

2.000 - 2.400 

2.194 

43.88 - 44.01 

43.94 

23.2 

(5.38 

± 

0.51) 

X 

10” 5 

2.000 - 2.400 

2.194 

44.01 - 44.13 

44.07 

30.3 

(5.15 

± 

0.43) 

X 

10 ” s 

2.000 - 2.400 

2.194 

44.13 - 44.26 

44.19 

33.3 

(3.63 

± 

0.30) 

X 

10 ” s 

2.000 - 2.400 

2.194 

44.26 - 44.38 

44.32 

25.3 

(2.33 

± 

0.22) 

X 

10” 5 

2.000 - 2.400 

2.194 

44.38 - 44.51 

44.44 

27.4 

(1.88 

± 

0.18) 

X 

10 ” s 

2.000 - 2.400 

2.194 

44.51 - 44.63 

44.57 

19.9 

(1.19 

± 

0.14) 

X 

10 ” s 

2.000 - 2.400 

2.194 

44.63 - 44.80 

44.72 

17.0 

(5.97 

± 

0.78) 

X 

10 -® 

2.000 - 2.400 

2.194 

44.80 - 44.96 

44.88 

9.0 

(2.44 

± 

0.44) 

X 

10 -® 

2.000 - 2.400 

2.194 

44.96 - 45.15 

45.06 

3.7 

(6.00 

± 

1.44) 

X 

10- 7 

2.000 - 2.400 

2.194 

45.15 - 45.39 

45.27 

2.8 

(2.94 

± 

1.09) 

X 

10- 7 

2.000 - 2.400 

2.194 

45.39 - 45.77 

45.58 

0.3 

(1.11 

± 

0.51) 

X 

10” 8 

2.000 - 2.400 

2.194 

45.77 - 46.02 

45.90 

0.9 

(1.22 

± 

1 .10) 

X 

10” 8 

2.400 - 3.000 

2.688 

42.71 - 43.21 

42.96 

3.6 

(5.25 

± 

1 .68) 

X 

10 ” s 

2.400 - 3.000 

2.688 

43.21 - 43.59 

43.40 

9.0 

(6.24 

± 

1 .12) 

X 

10 ” s 

2.400 - 3.000 

2.688 

43.59 - 43.72 

43.65 

8.5 

(8.04 

± 

1.61) 

X 

10 ” s 

2.400 - 3.000 

2.688 

43.72 - 43.84 

43.78 

7.4 

(4.45 

± 

0.89) 

X 

10” 5 

2.400 - 3.000 

2.688 

43.84 - 43.97 

43.90 

13.8 

(4.36 

± 

0.57) 

X 

10” 5 

2.400 - 3.000 

2.688 

43.97 - 44.09 

44.03 

14.0 

(3.09 

± 

0.37) 

X 

10 ” s 

2.400 - 3.000 

2.688 

44.09 - 44.22 

44.15 

22.2 

(3.07 

± 

0.30) 

X 

10 ” s 

2.400 - 3.000 

2.688 

44.22 - 44.34 

44.28 

19.3 

(1.98 

± 

0 .22) 

X 

10” 5 

2.400 - 3.000 

2.688 

44.34 - 44.47 

44.40 

23.9 

(1.61 

± 

0.16) 

X 

10” 5 

2.400 - 3.000 

2.688 

44.47 - 44.59 

44.53 

21.8 

(1.25 

± 

0.14) 

X 

10 ” s 

2.400 - 3.000 

2.688 

44.59 - 44.74 

44.67 

18.6 

(6.80 

± 

0.75) 

X 

10 -® 

2.400 - 3.000 

2.688 

44.74 - 44.87 

44.81 

12.3 

(4.16 

± 

0.62) 

X 

10 -® 

2.400 - 3.000 

2.688 

44.87 - 45.33 

45.10 

14.2 

(7.08 

± 

1.06) 

X 

10” 7 

2.400 - 3.000 

2.688 

45.33 - 45.58 

45.46 

0.6 

(3.99 

± 

1.36) 

X 

10” 8 

3.000 - 5.800 

4.215 

42.97 - 43.47 

43.22 

2.9 

(7.94 

± 

2.94) 

X 

10 -® 

3.000 - 5.800 

4.215 

43.47 - 43.75 

43.61 

3.7 

(7.03 

± 

1.76) 

X 

10 ”® 

3.000 - 5.800 

4.215 

43.75 - 43.97 

43.86 

8.5 

(9.67 

± 

1.74) 

X 

10 ”® 

3.000 - 5.800 

4.215 

43.97 - 44.10 

44.03 

5.9 

(5.56 

± 

1 .11) 

X 

10 ”® 

3.000 - 5.800 

4.215 

44.10 - 44.28 

44.19 

13.2 

(4.85 

± 

0 .68) 

X 

10 ”® 

3.000 - 5.800 

4.215 

44.28 - 44.41 

44.35 

5.5 

(1.70 

± 

0.31) 

X 

10 ”® 

3.000 - 5.800 

4.215 

44.41 - 44.53 

44.47 

8.3 

(1.87 

± 

0.30) 

X 

10 ”® 
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Figure [5] shows the variation of the smoothed 
two-power law parameters as a function of red- 
shift. The normalization is the XLF at log Lx = 
44 and we see the familiar feature of growth (with 
increasing redshift) between z = 0 and z ~ 1 , fol¬ 
lowed by a plateau, then a decline at z > 3. A 
notable feature is a sudden drop of the low lumi¬ 
nosity slope 71 from z < 0.6 (71 ~ 1 ), to z > 0.6 
(7 ~ 0.5). No significant change of slope is ob¬ 
served in each side of this redshift. An F-test com¬ 
parison between the best-fit models for the case 
where 71 is constant over redshift versus the case 
where 71 changes at z = 0.6 shows that the proba¬ 
bility that the former is accepted is only 7 x 10 -5 . 
The high luminosity slope is consistent with being 
constant at 72 ~ 2.7, except for two data points 
at z ~ 0.7 and z ~ 1.1, where 72 is significantly 
smaller, although it may be caused by an observa¬ 
tional bias as discussed in Sect. [3 


6.2. Evolution of the Space Density in Lu¬ 
minosity Classes 


In this section, we investigate the evolution 
of the AGN space density in different luminosity 
classes as a function of redshift. We have made the 
fit to an analytical expression at the center of each 
of the four luminosity classes: 42.0-43.0, 43.0- 
44.0, 44.0-45.0, and 45.0-46.0 in logL rmx . The 
XLF at the central luminosity of the bin logL c , 
which is defined as the mean of the minimum and 
maximum bounds of log L x bin can be expressed 
as: 


d${L c ,z) = AZ=0 
d log L x Lc 


e d (z,L c ), 


(24) 


where the normalization Af^° is the XLF value at 
the central luminosity of the class at z = 0 . 

For each luminosity class, we have made 
a maximum-likelihood fit to the the redshift- 
dependent evolution parameters in the evolution 
factor expressed as a three-segment power-law 
model Eq. [T51 evaluated at Lx = L c . 

During the fitting process, the luminosity de¬ 
pendence of the XLF within the shell is fixed to 
the best-fit LDDE case. In each class, there are 
six fitting parameters: A*j=°, pi, Zbi, P 2 , ~b 2 and 
P 3 . As in the case of redshift shell divided sam¬ 
ples, maximum-likelihood fits have been made for 
each luminosity class using Eq. [Tj where the nor¬ 
malization is also a fitting parameter. Thus the 



Fig. 5.— The parameters of the smoothed two 
power-law fits to individual redshift shells are plot¬ 
ted as a function of redshift. See labels. The error 
bars are ler corresponding to AL = 1. The up¬ 
ward arrow of in panel (d) corresponds to a 90% 
lower limit of 72 . 
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Table 4 —Continued 


z-range 

2c 

log Lx-range a 

logLx,c a 

•Wobj 

d& a ,b 

d log L x 

3.000 - 5.800 

4.215 

44.53 - 

■ 44.75 

44.64 

19.5 

(1.55 ±0.19) x 10 -6 

3.000 - 5.800 

4.215 

44.75 - 

■ 45.04 

44.90 

14.4 

(5.38 ±0.70) x 10 -7 

3.000 - 5.800 

4.215 

45.04 - 

■ 45.29 

45.17 

4.9 

(1.25 ±0.34) x lO" 7 


a Units- dlo g\ x : [h% Mpc 3 dex 1 ], L x : [h 70 2er S s 1 ]- 
b Errors show the 68% confidence range for one parameter (AC < 1). 


reported normalization errors contain the effects 
of correlations with other parameters. The result¬ 
ing best-fit parameters and errors are reported in 
Table 0 

One limitation of dividing the sample into lu¬ 
minosity classes is that we do not know the intrin¬ 
sic luminosity Lx of each object. We only know 
the observed flux and redshift. For the fit of the 
luminosity class (logLx,min,logLx,max), we select 
objects that fall into the log Lx,min < logLx.obs < 
log Lx.max range. Thus the fitting in each lumi¬ 
nosity class uses some objects that are outside of 
the log Lx, min < logL x ,obs < log Lx, max range, 
while others that fall into this range are not used. 
The fitting process using Eq. Improperly takes care 
of the expected number of AGNs in the observed 
luminosity range calculated based on the intrinsic 
luminosity XLF model. 

Figure 0 shows the evolution of the AGN num¬ 
ber density and emissivity in each class for the 
full treatment. The binned number densities 
and emissivities, calculated with the A^ 0 bs/AGidi 
method using the luminosity-class divided mod¬ 
els, are shown with 68% Poisson errors. The total 
emissivities of the 42 < logLx[/i7o ergs -1 } < 46 
range are also plotted as a function of redshift. 
The model number densities and emissivities are 
also shown for the best-fit LDDE model in Table [2] 
as well as class-divided fits. The evolution of total 
emissivity in 42 < log Lx < 46 is also plotted. Ta¬ 
ble [6] tabulates the binned Wibs/Anidi-estimated 
values evaluated at the central luminosities L c 
(instead of number densities integrated over the 
luminosity class, which are plotted in Fig. , as 
well as emissivities with best fit models. 

The results in Table [5] show that the first break 
redshift Zbi increases with luminosity with z\>\ ~ 
1.1 at log Lx ~ 44[h 7 o ergs -1 ]. The second break 


redshift does not vary significantly and is consis¬ 
tent with z b2 ~ 2.7 at log Lx > 43 [L 70 ergs -1 ]. 
The slope after the second break is consistent with 
p 3 —5.6. These have been fed back as fixed pa¬ 
rameters for a refinement of the global expression 
in Sect. 0 For the 42 < log Lx < 43 class, the 
values of Zb2 and p 3 could not be determined, be¬ 
cause the sample does not extend to the second 
break z^ 2 ~ 2.7. For the 45 < log Lx < 46 class, 
the low number of objects has limited our analy¬ 
sis. While we obtain best fit values of Zbi ~ 1-5 
and Zb2 ~ 3, their ler errors are essentially un¬ 
constrained. Thus we fixed the 2b 1 = 1.2 and 
Zb2 = 2.6 based on the results of the fit to the 
44 < log Lx < 45 class. We have also attempted 
a two-segment law fit (one break redshift) for this 
class, but the three segment form with the two 
break redshifts fixed at these values gives a better 
fit. 

While the value of the parameter pi increases 
with luminosity, the measurements are over differ¬ 
ent redshift ranges for different luminosity classes 
and should not be compared with one another. 
The p 2 parameters are also measured over some¬ 
what different redshift ranges. However, it is clear 
from Fig. [6] that the slope between Zbi and Zb2 
goes from negative to positive (as a function of 
log(l + 2)) with luminosity. The full binned XLF 
divided in luminosity classes is shown in Table 0 
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Table 5 

Best-fit evolution parameters for each luminosity class 


log Lx-range a 

log L c a 

Az=0a 

Lc 

Pi a 

z b i a 

P2 a 

*b2 a 

P3 a 

42.0 - 43.0 

42.5 

(8.3 ± 1.0) x 10" 5 

3.7 ±0.3 

0 65 - *" 0 04 
U.OO—003 

_o 4+ 0 - 5 



43.0 - 44.0 

43.5 

(6.0 ±0.4) x 10" 6 

4.7 ±0.2 

U.OO_Q 02 

-l-0la2 

2.69l®;“ 

-5.9+°-° 

44.0 - 45.0 

44.5 

(6.8 ±0.9) x lO" 8 

r O + 0.2 

°-°-0.3 

1 41 +0.09 

0-4lg;i 

2-63t“;(® 

-5.41H 

45.0 - 46.0 

45.5 

(8l®) x 10" 11 

7.3 ±0.8 

1.40* 

2.2 ±0.9 

2 .6* 

-4.7±i;l 


a The same units and error definitions as in Table [3] apply. 
''Fixed parameters are indicated by an asterisk (’*’). 




1+Z 1+Z 


Fig. 6.— (a) The space density of AGNs as a function of redshift in different luminosity classes as labeled 
in panel (b). Densities from the luminosity-class divided fits as well as those from the LDDE models are 
overplotted with solid lines and dotted lines respectively, (b) The 2-10 keV X ray emissivities are plotted 
instead of number densities. Also the total emissivity curves are added. Only models from the luminosity- 
divided fits are plotted in this panel. (Color versions of these figures are available in the electronic version.) 
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Table 6 

Full binned XLF divided in luminosity classes 


z -range 

Zc 

log L x - 

range a 

logL x ,c a 

Nobj 



d& a 

d Lx 

,b 


0.015 - 

0.089 

0.051 

42.00 - 

43.00 

42.50 

51.6 

(1.01 

± 

0.13) 

X 

10' 

0.089 - 

0.190 

0.138 

42.00 - 

43.00 

42.50 

21.2 

(1.02 

± 

0.13) 

X 

10' 

0.190 - 

0.276 

0.232 

42.00 - 

43.00 

42.50 

38.6 

(1.86 

± 

0.18) 

X 

10' 

0.276 - 

0.369 

0.322 

42.00 - 

43.00 

42.50 

51.4 

(2.38 

± 

0.20) 

X 

10' 

0.369 - 

0.468 

0.418 

42.00 - 

43.00 

42.50 

42.9 

(2.13 

± 

0.19) 

X 

10' 

0.468 - 

0.575 

0.521 

42.00 - 

43.00 

42.50 

57.3 

(3.29 

± 

0.25) 

X 

10' 

0.575 - 

0.689 

0.631 

42.00 - 

43.00 

42.50 

61.3 

(4.67 

± 

0.34) 

X 

10- 

0.689 - 

0.811 

0.749 

42.00 - 

43.00 

42.50 

45.6 

(3.60 

± 

0.29) 

X 

10“ 

0.811 - 

0.942 

0.875 

42.00 - 

43.00 

42.50 

29.2 

(2.70 

± 

0.26) 

X 

10- 

0.942 - 

1.083 

1.011 

42.00 - 

43.00 

42.50 

27.5 

(3.30 

± 

0.33) 

X 

io- 

1.083 - 

1.299 

1.188 

42.00 - 

43.00 

42.50 

16.9 

(1.86 

± 

0.26) 

X 

io- 

1.299 - 

1.615 

1.452 

42.00 - 

43.00 

42.50 

14.7 

(1.59 

± 

0.24) 

X 

io- 

1.615 - 

3.744 

2.522 

42.00 - 

43.00 

42.50 

12.0 

(4.18 

± 

0.67) 

X 

10“ 

0.015 - 

0.089 

0.051 

43.00 - 

44.00 

43.50 

165.3 

(6.56 

± 

0.51) 

X 

io- 

0.089 - 

0.168 

0.128 

43.00 - 

44.00 

43.50 

22.7 

(1.16 

± 

0.15) 

X 

io- 

0.168 - 

0.252 

0.209 

43.00 - 

44.00 

43.50 

30.3 

(1.69 

± 

0.19) 

X 

io- 

0.252 - 

0.343 

0.297 

43.00 - 

44.00 

43.50 

39.5 

(2.55 

± 

0.24) 

X 

IO" 

0.343 - 

0.440 

0.391 

43.00 - 

44.00 

43.50 

34.0 

(2.33 

± 

0.22) 

X 

io- 

0.440 - 

0.545 

0.492 

43.00 - 

44.00 

43.50 

56.4 

(3.50 

± 

0.26) 

X 

IO” 

0.545 - 

0.657 

0.600 

43.00 - 

44.00 

43.50 

75.7 

(4.42 

± 

0.28) 

X 

IO" 

0.657 - 

0.777 

0.716 

43.00 - 

44.00 

43.50 

140.4 

(8.33 

± 

0.37) 

X 

io- 

0.777 - 

0.906 

0.840 

43.00 - 

44.00 

43.50 

147.0 

(1.03 

± 

0.05) 

X 

IO" 

0.906 - 

1.044 

0.974 

43.00 - 

44.00 

43.50 

171.3 

(1.16 

± 

0.05) 

X 

IO" 

1.044 - 

1.192 

1.117 

43.00 - 

44.00 

43.50 

126.2 

(9.00 

± 

0.40) 

X 

io- 

1.192 - 

1.351 

1.270 

43.00 - 

44.00 

43.50 

119.2 

(9.54 

± 

0.44) 

X 

IO" 

1.351 - 

1.522 

1.435 

43.00 - 

44.00 

43.50 

99.1 

(9.13 

± 

0.43) 

X 

io- 

1.522 - 

1.704 

1.611 

43.00 - 

44.00 

43.50 

83.4 

(9.25 

± 

0.50) 

X 

10' 

1.704 - 

1.901 

1.801 

43.00 - 

44.00 

43.50 

49.9 

(6.77 

± 

0.43) 

X 

io- 

1.901 - 

2.111 

2.004 

43.00 - 

44.00 

43.50 

44.7 

(7.66 

± 

0.55) 

X 

io- 

2.111 - 

2.336 

2.222 

43.00 - 

44.00 

43.50 

32.2 

(6.79 

± 

0.59) 

X 

IO" 

2.336 - 

2.578 

2.455 

43.00 - 

44.00 

43.50 

27.7 

(7.25 

± 

0.70) 

X 

io- 

2.578 - 

2.838 

2.706 

43.00 - 

44.00 

43.50 

19.5 

(6.56 

± 

0.79) 

X 

IO" 

2.838 - 

5.309 

3.921 

43.00 - 

44.00 

43.50 

20.1 

(1.10 

± 

0.13) 

X 

IO" 

0.033 - 

0.108 

0.070 

44.00 - 

45.00 

44.50 

52.1 

(8.71 

± 

1.22) 

X 

IO" 

0.108 - 

0.282 

0.192 

44.00 - 

45.00 

44.50 

20.3 

(1.45 

± 

0.29) 

X 

10- 

0.282 - 

0.422 

0.350 

44.00 - 

45.00 

44.50 

21.9 

(4.74 

± 

0.76) 

X 

10- 

0.422 - 

0.590 

0.504 

44.00 - 

45.00 

44.50 

22.2 

(7.05 

± 

1.06) 

X 

10” 

0.590 - 

0.705 

0.646 

44.00 - 

45.00 

44.50 

22.4 

(1.55 

± 

0.23) 

X 

10” 

0.705 - 

0.829 

0.766 

44.00 - 

45.00 

44.50 

25.1 

(1.82 

± 

0.24) 

X 

10” 

0.829 - 

0.962 

0.894 

44.00 - 

45.00 

44.50 

43.4 

(3.20 

± 

0.29) 

X 

10- 

0.962 - 

1.104 

1.032 

44.00 - 

45.00 

44.50 

50.5 

(3.96 

± 

0.32) 

X 

10- 

1.104 - 

1.256 

1.179 

44.00 - 

45.00 

44.50 

71.3 

(5.95 

± 

0.40) 

X 

10- 

1.256 - 

1.420 

1.337 

44.00 - 

45.00 

44.50 

67.9 

(7.03 

± 

0.46) 

X 

10- 

1.420 - 

1.595 

1.506 

44.00 - 

45.00 

44.50 

107.7 

(1.33 

± 

0.07) 

X 

10- 

1.595 - 

1.784 

1.688 

44.00 - 

45.00 

44.50 

76.6 

(1.00 

± 

0.06) 

X 

10- 

1.784 - 

1.985 

1.883 

44.00 - 

45.00 

44.50 

101.5 

(1.45 

± 

0.07) 

X 

10” 

1.985 - 

2.202 

2.092 

44.00 - 

45.00 

44.50 

99.7 

(1.57 

± 

0.08) 

X 

10- 

2.202 - 

2.434 

2.316 

44.00 - 

45.00 

44.50 

78.2 

(1.40 

± 

0.08) 

X 

10- 

2.434 - 

2.683 

2.556 

44.00 - 

45.00 

44.50 

68.3 

(1.41 

± 

0.08) 

X 

10- 
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Figure [G^a) clearly shows a shift of the number 
density peak with luminosity up to log Lx ~ 44.5 
in the sense that more luminous AGNs (QSOs) 
peak earlier in the history of the universe, while 
the low luminosity ones arise later. Also, there 
is a clear decline of the derived space densities at 
all luminosity classes. A notable feature is a clear 
plateau structure at intermediate redshifts at in¬ 
termediate luminosities, which we have been able 
to trace in detail thanks to the addition of the 
COSMOS dataset. This has been reflected in the 
total emissivity curve, which stays almost constant 
between 1 < z < 2.7. The slope between the first 
break Zbi and the second break Zb 2 changes from 
negative to positive. 

7. Discussion 

7.1. Sensitivity to pz-PDZs and Spectral 
Models 

For comparison, we show the parameters for the 
case where we only take the best-fit single photo¬ 
metric redshift (PZ) per object, instead of full pz- 
PDFs, under “No PDZ” for the global expression 
(Table EJ . The fitted parameters are not signifi¬ 
cantly different between the two cases except £> 244 , 
which is the slope between the first break and the 
second break redshifts at log Lx = 44. The dif¬ 
ference in the normalizations are more than the 
errors given, but these errors are underestimates, 
since, for the global expression, we use Eq. [5] and 
correlations with other parameters have not been 
taken into account, as explained in Sect. 14.11 and 
15.21 To observe the effects more closely, we plot 
the ratio of the XLF estimates based on the full 
PDF and no PDF samples for each redshift shell 
in Fig. |T| The same are plotted for the num¬ 
ber density evolution in each luminosity class in 
Fig. [HI Upon the examining the difference be¬ 
tween the best-fit models in these figures, the ef¬ 
fects of the pz-PDF seem to be as large as ~ 30% 
at some locations (blue dashed lines). However, 
comparing the data points from the A r 0 b s /A r m di 
(black filled hexagonal versus blue open triangle 
data points in these figures), we see that there is 
no significant difference and the apparent discrep¬ 
ancy of the models are at locations that are not 
well-constrained by the data. Therefore the ap¬ 
parent deviation is not caused by the pz-PDFs, 
but by the different convergence of of the model 


parameters at the level of a fraction of statisti¬ 
cal errors. This is an optimistic case, because 
in the COSMOS survey, where the majority of 
the photometric redshifts with the pz-PDZ came 
from, has exceptionally high quality photometric 
data, achieving an ac curacy of Az/(1 + z ) « 0.015 
with ~ 5% outliers ( Salvato et all 2011 1. At the 
time of our analysis, full pz-PDFs were not avail¬ 
able for the CLANS, CLASXS, CDF-N and CDF- 
S data in our disposal, in which the spectroscopic 
redshift completeness ranges from 60-80%. These 
fractions are similar to those of the LH and COS¬ 
MOS dataset. For these, we take best-fit photo¬ 
metric redshift values, or for CDF-S sources, if 
the in formation about the second peak is avail¬ 
able, ( Luo et al. 20ld ; Xue et al. 2011 1. we give 
a weight of 0.5 in each peak for the full analy¬ 
sis. The photometric redshifts of the CDF-S have 
a normalized median absolute deviation ( NAM D) 


of onmad = 0.059 with outliers of 8 . 6 % (Luo et al 


2010ll . indicating a similar quality to those of COS¬ 


MOS. Those of the CLANS, CLASXS and CDF-N 
surveys use fewer photometric bands and the un¬ 
certainties are as la rge as a v J{\ + z) ~ 0.16, but 
with ~ 6 % outliers ( Trouille et al. 2008 1. The un¬ 
certainties of these PZ have not been taken into 
account in our full analysis. For a further test, 
we create fake PDZ’s for objects with photomet¬ 
ric redshifts (except those with a second peak) in 
these surveys with six delta functions. We assume 
a Gaussian profile centered at the cataloged PZ 
value z p for the main pz-PDF peak with an ap¬ 
propriate width corresponding to the onamd or 
<7 values for the each survey. In our fake PDZ, 
the main peak is represented by three delta func¬ 
tions, one at the center (representing the central 
1/2 of the PDF main peak), the others at the 
12.5% and 87.5% percentiles of the profile (repre¬ 
senting the upper and lower 1/4 of the PDF main 
peak). The central one has a weight that is half 
of 0.5(1 — outlierfraction) and the other two each 
has a weight that is half of the central one. Also 
added are three delta function to represent out¬ 
liers. Each of the three has a weight of 1/3 of the 
outlier fraction and the redshift is taken randomly 
from spectroscopic redshifts of X-ray sources that 
have fluxes within a factor of three of that the 
object. We keep using the real PDZs for the COS¬ 
MOS and LH samples. We make global fit the 
sample and all the fitted parameters are within a 
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Table 6 —Continued 


z-range 


log L x - 

range a 

log L XjC a 

Aobj 



a 

d L x 

l,b 



2.683 - 

2.950 

2.814 

44.00 - 

45.00 

44.50 

51.7 

(1.20 

± 

0.08) 

X 

10- 

-5 

2.950 - 

3.236 

3.091 

44.00 - 

45.00 

44.50 

29.2 

(7.85 

± 

0.73) 

X 

io- 

-6 

3.236 - 

5.749 

4.347 

44.00 - 

45.00 

44.50 

40.1 

(1.96 

± 

0.16) 

X 

io- 

-6 

5.749 - 

7.666 

6.648 

44.00 - 

45.00 

44.50 

0.0 

< 

1 

.1 x 10~ 

-6 


0.173 - 

0.470 

0.313 

45.00 - 

46.00 

45.50 

8.1 

(1.03 

± 

0.36) 

X 

io- 

-9 

0.470 - 

1.190 

0.794 

45.00 - 

46.00 

45.50 

6.9 

(4.18 

± 

1.21) 

X 

io- 

-9 

1.190 - 

1.480 

1.330 

45.00 - 

46.00 

45.50 

8.6 

(4.45 

± 

1.20) 

X 

10- 

-8 

1.480 - 

2.000 

1.728 

45.00 - 

46.00 

45.50 

10.5 

(5.67 

± 

1.13) 

X 

10- 

-8 

2.000 - 

2.600 

2.286 

45.00 - 

46.00 

45.50 

7.9 

(5.93 

± 

1.19) 

X 

10- 

-8 

2.600 - 

3.200 

2.888 

45.00 - 

46.00 

45.50 

7.7 

(8.42 

± 

1.77) 

X 

10- 

-8 

3.200 - 

5.310 

4.148 

45.00 - 

46.00 

45.50 

4.6 

(2.13 

± 

0.55) 

X 

10- 

-8 


a Units — dlo g\ x : [h 3 70 Mpc 3 dex x ], L x : [/i 70 2 ergs J ], 

b Errors show the 68% confidence range for one parameter (AC < 1). The symbol 
’<’ shows a 90% upper limit. 


fraction of ler errors from the original full analy¬ 
sis. Thus we conclude that the uncertainties of the 
PZ’s have minimal net effect to our XLF analysis. 


In the literature, it has been often assumed 
for the “hard” (2 < L[keV] < 7-10) sur¬ 
veys that absorption is negligible and the HXLF 
is calculated without taking absorption into ac¬ 
count, sim£ly_by_assuming a power-law spectra 
e.g. 


Silverman et al.l l2008t Yencho et al. 20091 


Aird et all 2010h ~ For comparison, we also com 


pute the case where all the AGNs have simple a 
r = 1.8 power-law. The results of this case, with 
PDFs of the photometric redshifts are taken into 
account, are listed under ‘T = 1.8” in Tableland 
the comparisons with the “Full” case are overplot¬ 
ted in Figs. [7] and [8] As seen in these figures, the 
XLFs based on the simple spectral assumption of 
r = 1.8 without absorption correction is subject 
to underestimates of the XLF at low luminosities 
(log Lx 44) at low redshifts (z < 1), by up to 
50%. If one uses a sample selected at E > 5 keV, 
however, this systematic error is expected to be 
reduced dramatically (S. Fotopoulou et al. 2015, 
submitted to A&A). 


In principle, the results depend on the IVh func¬ 
tion and its evolution. Using the 14-195 keV se¬ 
lected Swift BAT AGNs with detailed X-ray spec¬ 
troscopy at E < 10 keV and detailed modeling of 
selection effects, the Ah function derived by U14 
is robust, especially at low redshifts, and in agree¬ 


ment with that of Burlon et ah (2011) derived us¬ 
ing a larger sample of Swift BAT AGNs. The 
fraction of absorbed AGNs among CTN AGNs is 
obser ved to increase with redshift (e.g. Hasi nger 
20081 ). However, it is not clear whether the in¬ 
crease of absorbed fraction saturates at z ~ 2 or 
keeps increasing to higher redshifts, and Hasinger 
(120081) commented that a model where the ab¬ 
sorption fraction keeps increasing at 2 > 2 is also 
marginally acceptable. In our analysis above, we 
have used the U14 result, where the fraction of 
22 < log Ah < 24 AGNs among CTN AGNs at 
log L x = 43.75 grows as oc (1 + z) 0A5 at (z < 2.0) 
and stays constant at higher redshift. For a sensi¬ 
tivity check, we assume that the absorbed fraction 
grows as oc (1+z) 0 ' 32 up to our sample limit, where 
the index has been chosen to have the same ab¬ 
sorbed fraction at z = 3. The results of the global 
fit under this assumption have not changed signif¬ 
icantly, and all the parameters are consistent with 
our “Full” case within a fraction of ler. This can 
be understood because, at high redshift, the 0.5- 
10 keV band, which is usually used for the spec¬ 
tral analysis or hardness ratio analysis to measure 
Ah, corresponds to rest-frame E > 1.5 keV at 
z > 2. Thus it is already not sensitive to mea¬ 
suring absorption. For the same reason, from the 
2-10 keV flux, we measure rest frame higher en¬ 
ergy luminosity, which has not been affected by 
absorption. Since we want to derive the intrinsic 
XLF, the measured X-ray luminosity is insensitive 
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Fig. 7. The comparison of best fit models in each redshift shell and the ^V 0 bs/-^mdi estimates. The best 
fit two power-law models of the case of the no pz-PDFs and T = 1.8 are divided by those of the full case. 
The iV 0 bs/-/V m di XLFs are also divided by the full model. See labels in the lowest rightmost panel for the 
meanings of line styles, symbols, and colors (electronic version). The data points of the no PDF and P = 1.8 
cases are shifted horizontally by +0.04 and —0.04 respectively in log Lx for display purposes. (Color version 
of this figure is available in the electronic version.) 
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to the details of the IVh function at high redshifts. 

7.2. Detailed behavior of the XLF 

Thanks to the addition of the COSMOS 
dataset, we have been able to trace the detailed 
behavior of the XLF in the intermediate red- 
shift /intermediate luminosity range. We describe 
a number of notable features. 


• The low luminosity slope 71 flattens sud¬ 
denly at z ~ 0.6. The slopes stay al¬ 
most constant within each of the z < 0.6 
and z > 0.6 regimes, having 71 « 1 and 
~ 0.5 respectively. This sudden flattening of 
the slope practically excludes ILDE/LADE 
models, where the shape of the XLF stays 
unchanged over redshift, as a global expres¬ 
sion. The high redshift faint end slope is 
consistent with those of optically selected 
QSOs at 1 < z < 3.5 (71 ~ 0.45) mea¬ 
sured bv lBongiorno et all (20071) and z ~ 4 
QSOs bv llkeda et al.l (l201lh (71 = 0.7+0.2). 
Note, how ever, that the fa i nt end slope mea¬ 
sured by lOlikman et al. ( 2010 h is signifi¬ 
cantly steeper (71 = 1.3 ± 0.2 at z ~ 4). 


• The high luminosity slope stays almost con¬ 
stant (72 ~ 2.7), with the exception of two 
redshift shells at 0.6 < z < 0.8 and 1.0 < 
z < 1 . 2 , which have flatter slopes 72 « 2 . 0 . 
For the shell between these (0.8 < 2 < 1.0), 
72 is consistent with those for the rest of the 
shells. We note that the high luminosity end 
of the XLF is poorly sampled at these red- 
shifts and the apparently low 72 values may 
be caused by the higher weights at luminosi¬ 
ties just above the break and less weight at 
highest luminosities. Thus there is no evi¬ 
dence for systematic change of the high lu¬ 
minosity slope with redshift. The upcoming 
eROSITA mission will provide better sam¬ 
pling of this regime, which will enable us to 
see if the high end slope is constant over all 
redshifts. 


• The addition of the new data including 
COSMOS revealed new details about AGN 
downsizing. The number density/emissivity 
growth curves have two distinctive breaks, 
one at z^i ~ 1 and the other at 27,2 ~ 2.7, 
with a “plateau” between them. In any 



1+Z 


Fig. 8 .- The comparison of best fit models and 
the -/Vobs/Amdi estimates in each luminosity class. 
The best fit three segment power-law models of No 
pz-PDFs and T = 1.8 cases are divided by those 
of the full case. The N 0 t, s /N m< ii XLFs are also di¬ 
vided by the full model. The data points of the no 
pz-PDF and T = 1.8 cases are shifted horizontally 
by +0.007 and —0.007 respectively in log(l + z) 
for display purposes. See Fig. [7] for the meanings 
of line styles, symbols, and colors. (Color version 
of this figure is available in the electronic version.) 
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luminosity class, a rapid rise of number den¬ 
sity with cosmic time is observed at Zb2 (ex¬ 
cept the lowest luminosity bin) and a rapid 
decline is observed after Zbi- At low lu¬ 
minosities, (log Lx < 44), the peak is at 
Zbi and high luminosities, (log Lx < 44), 
Zb2 becomes the peak. Also the position of 
the first break moves from lower to higher 
redshifts with luminosity up from z ~ 0.6 to 
z ~ 1.4 up to log Lx < 44.5. No evidence for 
a change in Zb2 with luminosity is observed. 


Our results are in excellent agreement with 
those of U14 as seen in Fig. [9](a) . Small excesses 
at the level of < 0.2 dex at a few data points can 
be well explained by cosmic variance. Especially 
our afew x 10 15 < Sx.obs ^ 10 -14 ergs -1 cm -2 
range is dominated by the COSMOS data, which 
may be affecte d by large scale structures seen as 
redshift spike s ( Gilli et al . 200!): iBrusa et al . 20 lfl 


Civano et ah 20121) . Other possible differences at 
the same level may be explained in the treatment 
of the probability density distribution of the pho¬ 
tometric redshifts. Our results at the highest 
redshift shell z > 3 as well as the redshift evo¬ 
lution factor (p.s) are consistent with the recent 


results bv lVito et al.1 ( 20141 ). 


7.3. Comparison with Theoretical Models 


There are a number of attempts to repro¬ 
duce/understand the down-sizing behavior of 
AGN evolution within cosmological simulations 
and semi-analytical models. In order to under¬ 
stand the implication of our results to the physi¬ 
cal scenarios on accretion processes by comparing 
observations and theory, we should consider which 
features of the theoretical predictions are a conse¬ 
quence of the assumed physical picture and which 
features are a consequence of adjusting free pa¬ 
rameters to match the AGN luminosity functions 
in various redshifts. 


In their attempt to reprod uce the faint end 
AGN LF, iDeeraf et al.l ( 2010h consid ered a BH 


growt h by infall of surrounding gas via Bondi. H. 


(1952) accretion, which some authors call “radio 
mode” or “hot halo mode”, and merging with 
other BHs in their simulations including hydro- 
dynamic gas component. They compared their 
results with luminosity-class divided number den¬ 


sity evolution curves of soft (Hasinger et aL 2005 ) 


and hard X-ray ( Ueda et al. 20031) . They recog¬ 
nize that their model overproduces number den¬ 
sities of X-ray selected AGNs at z > 1. We 
compare our number de nsity evolution with two 
semi-analytical models of lMarulli et al.1 (120081 ) and 
Fanidakis et al. ( 20121) . where we have been able 


to convert their published results to the den¬ 
sity curves in our 2-10 keV luminosity classes in 
Fig. EKa). This figure shows that the overproduc¬ 
tion is also seen in these semi-analytical models 
(SAMs). 

Since these authors did not provide the num¬ 
ber density curves in the 2-10 keV range, we ex¬ 
plain the derivation of the luminosity-class di¬ 
vided 2-10 keV number densi t y cu rves from their 
published figures. For Marulli et al.l ( 20081 ) , we 
use their Fig. 7 and take the “best” bolomet- 
ric luminosity function model. The bolometric 
LF has been converted to 2-10 keV XLF model 


using the bolometric correction by iHopkins et al 


(120071) and integrated over our luminosity classes. 


For (F anidakis et al. (2012[), we use their soft (0.5- 
2 keV) X-ray “before obscuration (total)” curves 
provided in their Fig. 18 (b) instead of their “af¬ 
ter obscuration (visible)” to compare with our 2- 
10 keV density evolution curves. We convert their 
“before obscuration” to those for 2-10 keV assum¬ 
ing a photon index of T = 1.9, which is representa¬ 
tive of unobscured AGNs. Since they use slightly 
different cosmological parameters, we also convert 
the density curves to our adopted cosmology to 
overplot in Fig. Ha). 

In both cases, we see overpredictions of AGN 
number density at highest redshift at low lumi- 
nosity bins. As AGN model components, the 
Marulli et al. ( 2008i) model took into account 


the “radio” mode accretion and merger driven 
“quasar” mode. The lFanidakis et al.1 (|2012l) model 
took into account the “radio” and the “star- 
burst” modes, where the “starburst mode” in¬ 
cludes merger-driven “quasar mode” and secu¬ 
l ar (d isk instability) components. Marulli et al 


(|2008l ) considered a number of AGN light curve 
models and searched for the one that fits we ll 
with the bolometric LF of Hookins et al. ( 200 7). 
thus had freedo m in adjusting to o bserva tions. On 
the other hand, [Fanidaki s et al. ([20120 had four 
free parameters to fit to the LF. Except for the 
duty cycle, which is adjusted to the global normal¬ 
ization of the LF, the parameters are physically 
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motivated. The Fanidakis et al. ( 2012h model 


overpredicts the number densities at high red- 
shift at all luminosity classes. In their compari¬ 
son with the 0.5-2 keV XLF, they took into ac¬ 
count the effects of luminosity dependence and 
redshift evolution of obscu red/type 2 AGN frac¬ 
tion f- 2 by Hasingeil (l2008h. For the redshift evo¬ 
lution of / 2 , Hasingeil ~ (j2008} ) proposed a model 
with / 2 oc (1 + z ) 0 62 up to ^ = 2.06 and which 
stays constant at hig h er re dshifts as preferred, 
while Fanidakis et al. ( 2012h used a marginally 
acceptable model with / 2 oc (1 + z) 0 ' 48 without 
saturation, and thus their obscured fraction keep 
increasing to higher redshifts. The latter form 
of the obscured fraction evolut io n wa s needed 
in the model of iFanidakis et aT ( 2012 ) for sat¬ 
isfactory consistencies with the optical QSO and 
soft XLF to their model. On the other hand, 
the high-z overprediction of the 2-10 keV number 
densities show that they indeed over predicted 
high redshift AGN population and over-corrected 
for the effects of obscuration to match with the 
0.5-2 keV number density/optical QSOs. The 
tendency of overprediction of high redshift low lu¬ 
minosity/mass population is also a problem with 
predicting stellar mass function of galaxies (e.g. 

IGuo et ali 1 2011 ). the reason for which is still un¬ 
clear. On the other hand, if their basic prediction 
on high redshift accretion is correct, the discrep¬ 
ancy may be an indication of a large population 
of highly Compton-thick AGNs at high redshifts. 

The high redshift over-prediction problem is al so 
present in the recent model bylSiiacki et al. (2014) 
at all luminosities. 

Another interesting model is bv lDraper &; Ballantvnel 
([2012), where they considered merger and secular 
modes of AGNs triggering. Their work was not 
directly based on cosmological simulations. In¬ 
stead they took simple formulae for the evolution 
of galaxy number density with stellar mass above 
some threshold value, gas fraction, and merger 
rate from the literature. They then considered sec¬ 
ular and merger driven AGN triggering followed 
by a parameterized model of luminosity decay. 

Their model 2-10 keV number density evolutions 
are overplotted with our data in Fig. [U)b), which 
has been derived from their Fig. 7. Since their Fig. 

7 shows curves for log Lx > 42, > 43, and > 44, 
we draw their model curves of log Lx = 42 — 43 
and 43 — 44 curves by subtracting their > 43 curve 


from > 42 one and > 44 curve from > 43 one re¬ 
spectively. For the 44—45 curve, we use their > 44 
curve, since the contribution of the log Lx > 45 
AGNs to the number density is only ~ 0.2% of 
> 44 AGNs for y 2 = 2.7. We see that their “sum” 
of the contributions of the two processes is in 
good agreement with our data. Since they have 
three free parameters to adjust in their AGN light 
curves, which are the timescale, slope of the de¬ 
cay, and the peak Eddington ratio, they had much 
freedom in adjusting to observed XLFs. Thus it is 
no surprise that their theoretical curves agree well 
with our observation. However, they had difficul¬ 
ties in reproducing the hard X-ray number density 
curve with only one mode, i.e., merger only or sec¬ 
ular only. Their summed curves show the flat-top 
(two break) structures, where the slope in the 
Zbi < z < Zb2 range increases from negative to 
positive with luminosity. 

8. Conclusions 

We have investigated the detailed behavior of 
the intrinsic 2-10 keV X-ray luminosity function 
(XLF) of active galactic nuclei (AGNs) from a 
combination of samples in various depths/width. 
We summarize our conclusions. 

1 . A total of ~ 3200 unique X-ray AGNs from 
a combination of wide to deep samples cov¬ 
ering the redshifts of 0.015 < z < 5.8 and 
six orders of magnitude in flux have been 
used to investigate the behavior of the XLF. 
In particular, we present here the XLF from 
X-ray sources from both of the XMM and 
Chandra COSMOS surveys, which comprise 
about 43% of our global sample. 

2. Recent Ah functions and realistic spec¬ 
tral templates are fully incorporated in the 
maximum-likelihood fitting to analytic ex¬ 
pressions as well as the A 0 b s /A m di-estimated 
binned intrinsic XLF. The photometric red¬ 
shift probability density distributions (pz- 
PDFs) of the COSMOS and Lockman Hole 
samples as well as secondary photometric 
redshifts of CDF-S are also considered. 

3. We present a full parametric luminosity- 
dependent density evolution (LDDE) model, 
and redshift-shell and luminosity class sep- 
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Fig. 9. — (a) Our luminosity-class divided num¬ 
ber density evolution curves (binned data points 
with errorbars from Fig. (6]i) are compared with 
the LDDE model fr om recent comparable work 
bv lUeda et al. ( 20141 ) (dotted lines). Also results 
from recent semi-analytca l models from literature 
from Marulli et al. (l2008h (dot-dashed lines) and 
Fanidakis et al. ( 2012h (dashed lines) are overplot¬ 


ted. Luminosity classes for the data points and 
model curves are indicated by colors as labeled. 
See text for detailed conversions of their published 
results to the 2-10 k eV density curves. Model by 


Marulli et al. ( 2008 1 does not reach the highest 


luminosity bin and therefore we are not able to 
plot log Lx = 45 — 46 bin. (b) The same num¬ 
ber d ensity curves are c o mpar ed with the model 
bv iDraper & Ballantvne ( 2012i l (solid lines) with 
separate contributions from merger (dotted lines) 
and secular processess (dashed lines). (Color ver¬ 
sions of these figures are available in the electronic 


version. 


arated models, as well as the binned XLF 
using the -/Vobs/Amdi estimator. 

4. The low luminosity end slope of the XLF 
flattens suddenly at z > 0.6. This behavior 
practically excludes any global expression of 
XLF evolution that assumes an unchanged 
XLF shape over cosmic time. 

5. We investigate the net effects of ignoring 
probability density distribution (pz-PDF) of 
photometric redshifts. The effect of includ¬ 
ing pz-PDF into the analysis does not alter 
the final XLF results significantly. 

6 . We investigate the effects of assuming a sim¬ 
ple AGN spectrum of T = 1.8. Under this 
simple assumption, the XLF is subject to un¬ 
derestimates by up to ~ 50% in the low red- 
shift low luminosity regime. 

7. The detailed behaviors of AGN down-sizing 
have been revealed. A clear two-break struc¬ 
tures has been revealed in the number den¬ 
sity evolution curves at every luminosity 
class above log Lx > 43. This behavior 
is qualitatively consistent with a two-mode 
AGN evolution involving major merger and 
secular processes. 

8 . Most current semi-analytical models of 
AGNs over-produce AGN number densities 
in the high redshift, low luminosity regime. 
Their semi-analytical treatment of accre¬ 
tion processes may have to be revised to 
reproduce our results. Alternatively, if their 
accretion scenario is correct, this may sug¬ 
gest a population of heavily Compton-thick 
AGNs in this regime. 
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